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ABSTRACT 

We use a sample of 90 spectroscopically-confirmed Lyman Break Galaxies with Ha measurements 
and SpitzerMIPS 24 fim observations to constrain the relationship between rest-frame 8 fxm luminosity 
(Lg) and star formation rate (SFR) for L* galaxies at z ~ 2. We find a tight correlation with 0.24 dcx 
scatter between Lg and Ha luminosity /SFR for z ~ 2 galaxies with 10 10 < LrR < 10 12 L Q . Employing 
this relationship with a larger sample of 392 galaxies with spectroscopic redshifts, we find that the UV 
slope /3 can be used to recover the dust attenuation of the vast majority of moderately luminous L* 
galaxies at z ~ 2 to within a 0.4 dex scatter using the local correlation. Separately, young galaxies with 
ages < 100 Myr appear to be less dusty than their UV slopes would imply based on the local trend 
and may follow an extinction curve that is steeper than what is typically assumed. Consequently, 
very young galaxies at high redshift may be significantly less dusty than inferred previously. Our 
results provide the first direct evidence, independent of the UV slope, for a correlation between UV 
and bolometric luminosity (Lboi) a t high redshift, in the sense the UV-faint galaxies are also on 
average less infrared and less bolometrically-luminous than their UV-bright counterparts. The L^ \- 
Luv relation indicates that as the SFR increases, Luv turns over (or "saturates") around the value of 
L* at z ~ 2, implying that dust obscuration may be largely responsible for modulating the bright-end 
of the UV luminosity function. Finally, dust attenuation is found to correlate with oxygen abundance 
at z ~ 2. Accounting for systematic differences in local and high-redshift metallicity determinations, 
we find that L* galaxies at z ~ 2, while at least an order of magnitude more bolometrically-luminous, 
exhibit ratios of metals-to-dust that are similar to those of local starbursts. This result is expected 
if high-redshift galaxies are forming their stars in a less metal-rich environment compared to local 
galaxies of the same luminosity, thus naturally leading to a redshift evolution in both the bolometric 
luminosity - metallicity and bolometric luminosity - obscuration relations. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: high-redshift — infrared: 
galaxies — ISM: dust, extinction 



1. INTRODUCTION 

Quantifying the total energetics and the time evolution 
of the dust properties of high-redshift galaxies requires 
an understanding of the extent to which the locally cali- 
brated relations between mid-infrared (rest-frame 8 /mo) 
and bolometric luminosity apply at high redshift. The 
Spitzer Space Telescope allows for a direct measure of 
the mid-IR dust emission from typical (L*) galaxies at 
z > 1.5 and has motivated a number of investigations 
of the correlation between mid-IR and bolometric lu- 
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minosity. The mid-IR emission from 5 — 8.5 /j,m (rest- 
frame) arises from the stochastic UV photo-heating of 
small dust grains and polycyclic aromatic hydrocarbons 
(PAHs; e.g., Puget & Leger 1989; Tielens et al. 1999). 
As such, this emission is found to correlate with the UV 
radiation from OB stars and hence the global star for- 
mation rate in nearby galaxies (e.g., Forster Schreiber 
et al. 2004; Roussel et al. 2001), albeit with some varia- 
tion depending on metallicity and ionizing hardness (e.g., 
Engelbracht et al. 2005; Hogg et al. 2005; Helou et al. 
2001; Normand et al. 1995. More recently, Kennicutt 
et al. (2009) demonstrate for the local SINGSs sam- 
ple of galaxies the general agreement between Balmer- 
decrement dust-corrected Ha star formation rates and 
those derived by combining tracers of obscured star for- 
mation, including IR and 8 /xm luminosity, and tracers 
of unobscured star formation (observed Ha luminosity). 

At high redshift, several studies suggest that on av- 
erage the ratios of mid to total infrared luminosity for 
luminous 24 /im-selected star-forming galaxies at z ~ 2 
are larger than those found for local galaxies with simi- 
lar bolometric luminosities (Papovich et al. 2007; Rigby 
et al. 2008). On the other hand, X-ray and radio stack- 
ing analyses suggest that the correlation between mid-IR 
and infrared luminosity for the less luminous but more 
typical galaxies at z ~ 2 (i.e., those with luminosities 
comparable to L^ v or L£ ol at z ~ 2; Reddy et al. 2008; 
Reddy & Steidel 2009) is consistent with the local rela- 
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tions (Reddy et al. 2006b; Reddy & Steidcl 2004). Taken 
together, these results imply a luminosity-dependence in 
the correlation between mid-IR and bolometric luminos- 
ity. 

While progress in connecting the dust emission of high- 
rcdshift galaxies to their bolometric luminosities can be 
informative, most z > 3 galaxies are selected via the effi- 
cient UV-dropout technique (Steidel et al. 1 995) and tend 
to be too faint and at too high redshift to be detected di- 
rectly via their dust emission. Thus, it has become com- 
mon to rely on UV-based inferences of the dust attenua- 
tion and bolometric star formation rates of z > 3 galaxies 
that require some assumption of how the rest-frame UV 
slope varies with extinction. The highest redshift where 
the correlation between UV-slope and dust extinction 
has been tested for large numbers of spectroscopically- 
confirmed galaxies is at z ~ 2. The results suggest that 
the local relation between UV-slope and extinction re- 
mains valid at z <~ 2 for galaxies with moderate bolo- 
metric luminosities similar to those of luminous infrared 
galaxies, or LIRGs (Reddy et al. 2006b; hereafter R06). 

Access to a third tracer of star formation, independent 
of the UV and IR emission, is required to establish a 
more secure footing for (1) the scaling between mid-IR 
emission and star formation rate and (2) the dependence 
of UV-slope on dust extinction. The most direct probe 
of young massive star formation is from HII recombina- 
tion nebular emission. In particular, the Ha line has been 
used traditionally as a star formation rate indicator given 
its accessibility in the optical window and the fact that it 
traces star formation on very short timescales (~ 10 Myr) 
and is less sensitive to extinction than UV emission (e.g., 
Kennicutt 1998; Brinchmann et al. 2004). Fortunately, 
Ha is still accessible with ground-based near-IR spectro- 
graphs for galaxies at z ~ 2. Using a sample of 114 
UV-selected galaxies at z ~ 2 with Ha spectroscopy, 
Erb et al. (2006b) demonstrate that the Ha inferred 
star formation rates agree well with those obtained from 
the UV after correcting both the Ha line and the UV 
continuum magnitude for extinction, and assuming that 
E(B - V) st eiiar ~ E(B - U) no buiar- In this paper, we use 
this Ha spectroscopic sample as the basis for quantify- 
ing independently the scaling between mid-IR emission 
and star formation rate for moderately-luminous galax- 
ies at z ~ 2. This scaling relation is then used to infer 
bolometric luminosities and dust obscuration and deter- 
mine how the UV slope varies with extinction for typical 
high-rcdshift galaxies. 

We begin by discussing in § 2 the rest-UV selection, fol- 
lowup optical and near-IR spectroscopy, and the MIPS 
observations, data reduction and photometry. The stel- 
lar population modeling of galaxies in our final samples 
is discussed in § 3. We then present the correlations 
between Ha, X-ray, and 8 /urn luminosity (and the con- 
nection between these quantities and the star formation 
rate) and comparison with other local and high-redshift 
relations (§4). The variation of UV slope with dust at- 
tenuation is addressed in § 5. The correlation between 
bolometric luminosity and dust attenuation and the im- 
plications of this relationship for the observed UV lumi- 
nosity are presented in § 6. Finally, in § 7 we present 
the results of a comparison between extinction and gas- 
phase metallicities of star-forming galaxies at z ~ 2 and 
demonstrate a close relationship between the two. We 



assume a Chabrier (2003) IMF unless stated otherwise 
and adopt a cosmology with Hq = 70 km s" 1 Mpc -1 , 
n A = 0.7, and ^ m = 0.3. 

2. SAMPLE 

2.1. Optical and Ha Observations and Spectroscopy 

The photometry and spectroscopic followup for the Ly- 
man Break Galaxy (LBG) survey at redshifts 1.4 < z < 
3.4 are described in Steidel et al. (2003, 2004); Adelberger 
et al. (2004). Briefly, galaxies were photometrically- 
selected using either the "BM," "BX," or LBG rest-UV 
color criteria. A subset of 1Z < 25.5 candidates that 
were spectroscopically-confirmed with Keck/LRIS obser- 
vations (Steidel et al. 2004) were targeted with NIRSPEC 
H or K-band spectroscopy to trace Ha (Erb et al. 2006c). 
In addition to the longslit Ha spectroscopy, we targeted 
the z = 2.300 redshift overdensity in the HS1700+64 
("Q1700") field (Steidel et al. 2005) with narrowband 
Br7 (2.17 /zm/0.04 /im) observations that trace Ha at 
that redshift. These images, along with JT s -band ob- 
servations, were obtained with the Wide-Field Infrared 
Camera (WIRC) on the Palomar 5 m telescope. De- 
tails on the Br7 imaging reductions, continuum subtrac- 
tion, and Ha flux determinations will be provided in a 
forthcoming paper (Steidel et al., in prep.). Briefly, we 
calibrated the Br7 image by adjusting its zeropoint un- 
til the difference between K s and Br7 magnitude was 
zero for objects of similar faintness in K s . Colors were 
measured in matched isophotal apertures with detection 
done at -fT s -band. We then measured the narrowband 
fluxes and subtracted the continuum based on the Br7- 
K s colors and knowledge of the Br7 filter bandwidth and 
shape. Photometric and zeropoint uncertainties result in 
an ~ 20% error in the derived line fluxes. 

This procedure resulted in 21 objects with > 5 a flux 
limits in Ha, ranging in flux from <~ 5 x 10~ 17 to <~ 
3 x 10~ 16 ergs s" 1 . Eleven of these 21 objects also have 
longslit Ha observations. Combining the longslit and 
narrowband Ha samples yields a total number of 116 
galaxies in 6 fields that have 24 /im coverage and are 
spectroscopically confirmed to lie at 1.5 < z < 2.6, where 
the 24 ^m data are sensitive to the rest-frame 8 ^m flux 
(§ 2.3). 

2.2. Correction for Slit Loss 

The primary uncertainty in the spectroscopically- 
derived Ha fluxes is slit loss. The overlap of 11 objects 
between the Ha longslit and narrowband samples in the 
Q1700 field allows us to quantify this source of error. To 
assess the degree of [Nil] contamination of the narrow- 
band fluxes, we constructed a composite spectrum of the 
11 objects with both longslit and narrowband observa- 
tions. From this composite, we determine a mean ratio 
of the [NII]-to-Ha flux of w 0.09 for these 11 objects. At 
the redshift of these galaxies (z — 2.300), the Br7 trans- 
mission curve will suppress the [Nil] line by a factor of 
0.88 relative to Ha. Therefore, we expect the average 
[Nil] line contamination of the narrowband flux to be 
~ 8%, translating to a 0.03 dex correction in log flux. 
After taking into account this 8% correction, the ratio of 
the narrowband Ha flux to the spectroscopic Ha flux for 
these 11 objects varies from 1.14 to 4.21, with a median 
value of 1.87 and a 1 a dispersion around this median 



3 



value of ~ 0.8 (i.e., a scatter of w 40%). This median 
value is similar to the average slit loss correction factor 
of 2.0 suggested by Erb et al. (2006b). For galaxies with 
near-IR spectroscopic observations only, we correct the 
Ha flux for slit losses by multiplying by the median value 
found above, 1.87. 

The sample of 11 galaxies with both longslit and nar- 
rowband observations have an Ha flux distribution that 
is undistinguished from the distribution for the 90 galax- 
ies that are used to measure the relation between 8 mi- 
cron luminosity (L& — i/L w (8jj l m)) and L(Ha). Over this 
flux range, we find no dependence of the slit loss correc- 
tion factor with the spectroscopic flux. The absence of 
any spectroscopic flux-dependent systematic in the slit 
loss correction, combined with the expectation that the 
amount of Ha flux lost in the spectroscopic observations 
depends sensitively on the relative placement and posi- 
tion angle of the slit on the galaxy, leads us to the con- 
clusion that the slit loss correction must be dominated 
by random error. 

2.3. MIPS Observations 

There are 6 fields in the LBG survey that have Spitzer 
MIPS 24 /im coverage, including the publicly available 
data in the GOODS-N and Westphal fields (Dickinson 
et al., in prep.; Chary et al., in prep.). We obtained 
observations in 4 additional LBG fields: Q1623, Q1700, 
Q2343, and Q1549 (see Reddy & Steidel 2009). Obser- 
vations in these 4 fields consisted of 30 sec frame times 
with 2 mapping cycles, where each cycle consists of a 
3x3 dither pattern with a step size of 1/8 the size of the 
MIPS array. The total exposure time in 3 of the 4 fields 
is « 11.4 hours, sufficient to detect ~ 15 /iJy sources at 
> 3 a (in the absence of confusion) within the inner 4' of 
each field. This is comparable to the depth of the 24 /<m 
data in the GOODS-N and Westphal fields. The Q1700 
MIPS data consist of « 20 hrs of observations and reach 
a 3 (J depth of « 10 /iJy in the absence of confusion. 

The 24 /im images are processed using a procedure 
similar to that described in R06. Briefly, the data are 
flatficlded and then combined using MOPEX (Makovoz 
& Marleau 2005). Photometry is extracted using PSF 
photometry with priors determined by the locations of 
objects detected in the IRAC images (IRAC data exist 
in all 6 fields). There are three sources of uncertainty in 
the 24 /mi fluxes. The first is background noise, which we 
determine by fitting PSFs simultaneously to random po- 
sitions around the target of interest and any other known 
sources. The dispersion in measurements of the back- 
ground level obtained in this way is similar to the disper- 
sion we would estimate by comparing the simulated and 
recovered fluxes of artificial objects added to the 24 /im 
images. The second source of uncertainty is Poisson error 
and is generally negligible compared to the background 
dispersion. The third source of uncertainty is confusion, 
which is reflected by the degeneracy between the PSF 
fit to the object of interest and the fits to any nearby 
sources. The uncertainty in the fit is computed simply 
from the covariance matrix for the N number of objects 
that are being fit simultaneously. The total (normalized) 
covariance in the fit (cr cov ) is a number between and 1, 
where the former indicates no covariance with nearby ob- 
jects and the latter indicates total covariance (e.g., such 
as the covariance between the fit of an object and itself). 



For most galaxies in the sample of 116, cr cov << 1. There 
are 17 out of 116 galaxies that have cr cov > 0.3. We con- 
sider these sources to be confused and remove them from 
our sample: this results in 99 galaxies with secure PSF 
fits. 

2.4. Identification of AGN 

Because we are interested primarily in the star-forming 
galaxies, we have removed AGN from the sample based 
on (1) the presence of prominent emission lines includ- 
ing Lya, CIV, and NV, or (2) an SED that behaves 
as a power-law between rest-frame ~ 2 /jm and 8 /im 
(§ 3). Among the 99 unconfused galaxies in our sam- 
ple, 9 are identified as AGN. Roughly half of these AGN 
also exhibited high ratios of [Nil] /Ha > 0.5. The frac- 
tion of AGN (« 9%) is somewhat larger than the ~ 5% 
found among typical LBGs to TZ = 25.5 (Reddy et al. 
2006a); the Ha spectroscopic sample is biased towards 
more K-bright galaxies relative to the typical LBG (Erb 
et al. 2006b) , and the frequency of AGN increases rapidly 
with increasing if-band brightness at z <~ 2 (Reddy et al. 
2005). 

2.5. Final Samples 

The sample used to measure the relationship between 
mid-IR luminosity and star formation rate consists of 90 
galaxies in the Ha sample (including spectroscopic and 
narrowband selected objects) that satisfy the following 
criteria: (1) 1.5 < ^spec < 2.6, (2) cr cov < 0.3, and (3) do 
not show AGN signatures based on optical emission lines 
or a power-law SED. Of these 90 galaxies, 29 are unde- 
tected at 24 /im to 3 a. We use a larger sample of 158 
UV-selected galaxies in the GOODS-N field that satisfy 
the three criteria specified above and have correspond- 
ing X-ray data that allow us to perform an X-ray stack- 
ing analysis (§ 4.2). Finally, the relationship between 
dustiness and UV spectral slope is investigated with a 
sample of 392 galaxies (from 5 of the 6 LBG fields with 
MIPS data) that satisfy the three criteria above and have 
had their photometry modeled using a stellar population 
analysis (§3). 

3. STELLAR POPULATION MODELING 

To provide a broader context for interpreting the ex- 
tinction properties, we have derived ages and stellar 
masses by modeling the broadband photometry of galax- 
ies in our sample. All of the galaxies with MIPS data 
have been imaged with IRAC, enabling us to use the 
full rest-frame UV through near-IR photometry to fit for 
their stellar populations. Previous efforts in this regard 
are described in Shapley et al. (2005), Erb et al. (2006c), 
and Reddy et al. (2006a). Here, we update these re- 
sults using the latest models of Chariot & Bruzual (in 
prep; CB08) that include the Marigo & Girardi (2007) 
prescription for the thermally-pulsating AGB (TP-AGB) 
evolution of low- and intermediate-mass stars. The ef- 
fect of this evolution on the inferred ages and stellar 
masses is then determined by comparing with our previ- 
ous results that employed the Bruzual & Chariot (2003) 
(BC03) models. 

3.1. Modeling Procedure 

For each galaxy, we considered a constant star- 
formation (CSF) model and exponentially-declining star 
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Fig. 1. — Comparison of the best- fit ages {left) and stellar masses {right) of 95 galaxies with Ha measurements before and after correcting 
the i^s-band flux for Ha line contamination. Points are color-coded according to observed Ha equivalent width. 



formation histories with characteristic time-scales r = 
10, 20, 50, 100, 200, 500, 1000, 2000, and 5000 Myr. We 
further considered a range of ages spaced roughly loga- 
rithmically between 1 and 5000 Myr, excluding ages older 
than the age of the universe at the redshift of each galaxy. 
In addition, allowed ages were restricted to be longer 
than the dynamical time-scale of ~ 70 Myr, as inferred 
from velocity dispersion and size measurements of z ~ 2 
LBGs (Erb et al. 2006b; Law et al. 2007). Reddening is 
taken into account by employing the Calzetti extinction 
curve and allowing E(B — V) to range between 0.0—0.6. 
In S 5.3 we also consider the effects of adopting the SMC 
extinction curve on the model results for the youngest 
galaxies in our sample. The model SED at each r and 
age combination is reddened, redshifted, and attenuated 
blueward of rest-frame 1216 A for the opacity of the 
IGM using the Madau (1995) prescription. The best-fit 
normalization of this model is determined by minimizing 
its x 2 with respect to the observed U n GlZ + J-K S +IRAC 
(3.6—8.0 /zm) photometry. This normalization then de- 
termines the star formation rate and stellar mass. The 
model (and normalization) that gives the lowest x 2 is 
taken to be the best-fit SED. We note the inherent de- 
generacy in this type of modeling. Typically there are 
several best-fit models that may adequately describe the 
observed photometry, even when the redshift is fixed to 
the spectroscopic value, though there is generally less 
variation in stellar mass than in the other parameters 
(r, age, E(B — V)) among these best-fit models (Shapley 
et al. 2005, 2001; Papovich et al. 2001; Sawicki & Yee 
1998). 

In the subsequent analysis, we have adopted the best- 
fit parameters obtained with a CSF model for several 
reasons. This model generally yields x 2 values similar to 
those obtained when t is allowed to vary freely. Addi- 
tionally, some of the more extreme star formation histo- 
ries are ruled out because the youngest galaxies in our 
sample cannot realistically have ages much less than the 
dynamical time-scale of ~ 70 Myr. Such extreme mod- 
els are also unlikely based on the presence of O star and 
Wolf-Rayet stellar features in the composite UV spectra 



of z ~ 2 galaxies, irrespective of age (e.g., Shapley et al. 
2005). For simplicity, we do not consider more complex 
multi-component or stochastic star formation histories 
as the data would not allow for discrimination between 
them as compared to simpler star formation histories. 
Finally, we do not use the SFRs from the SED fitting in 
the subsequent analysis, opting instead to use the direct 
tracers examined here, including UV, Ha, X-ray, and 
8 /jm luminosity. 

3.2. Effect of Ha Emission on the SEDs 

The parameters of most interest are the age and stel- 
lar mass, both of which are determined primarily by the 
strength of the Balmer and 4000 A breaks and the stel- 
lar continuum flux as traced by the near-IR and IRAC 
photometry at z ~ 2. Because the Ha emission line 
falls in the -K s -band at z ~ 2.0 — 2.6, it is prudent to 
determine what effect, if any, the line has on biasing 
our estimates of the ages and stellar masses. Figure 1 
summarizes the comparison of the best-fit ages and stel- 
lar masses (assuming constant star formation) obtained 
with and without correcting for the Ha contribution to 
the .Kg-band flux for 95 galaxies with Ha measurements. 
Roughly two-thirds of the sample of 95 galaxies show no 
change in best-fit age or stellar mass after correcting the 
-ftf s -band flux for Ha line contamination, and not surpris- 
ingly these galaxies generally have small Ha equivalent 
widths (Wobs ;$ 300 A). Even for those galaxies with 
larger Ha equivalent widths, the mean difference in age 
and stellar mass is relatively small compared to the un- 
certainty in these parameters as a result of the significant 
degeneracies inherent in SED fitting. 

The small differences in the ages and stellar masses be- 
fore and after correcting for Ha emission partly reflect 
the fact that these parameters are not solely constrained 
by the near-IR data: 89 of the 95 galaxies also have 
higher S/N IRAC data that are unaffected by line con- 
tamination at these redshifts. Also note that the age 
distribution for galaxies with the largest Ha equivalent 
widths (W ohs > 1000 A), including many of those se- 
lected from the narrowband data of the Q1700 field (§ 2), 
is not significantly different from that of galaxies with 
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Fig. 2. — Ratio of BC03 to CB08 derived stellar masses as a 
function of BC03 derived stellar mass for 1156 UV-selected galaxies 
with spectroscopic redshifts 1.40 < z < 3.50. Fractions indicate the 
percentage of objects in each age range that have CB08 masses that 
deviate more than 20% from the BC03 mass. 

smaller equivalent widths. The lack of a significant trend 
between Ha equivalent width and age for this sample has 
been noted by Erb et al. (2006b), and may be due to 
the stochasticity of the instantaneous SFR. In any case, 
given that the age and SFR distribution of galaxies with 
Ha measurements span roughly the entire range probed 
by UV-selected galaxies in general (Erb et al. 2006c), 
and that the correction for the Ha flux does not signif- 
icantly alter the best-fit ages and stellar masses for the 
majority of the sample, we are confident in our ability 
to constrain the SED fits of galaxies at z ~ 2.0 — 2.6 
even in the absence of direct measurements of the Ha 
line contamination in the AVband. 

3.3. Comparison between CB08 and BC03 

Figure 2 compares the BC03 and CB08 derived stellar 
masses for 1156 UV-selected galaxies with spectroscopic 
redshifts 1.40 < z < 3.50, color coded by the BC03- 
derived stellar age. The TP-AGB phase will become im- 
portant ~ 200 Myr after the onset of star formation. 
Anywhere from 50 — 60% of galaxies older than this have 
CB08 derived stellar masses that deviate by more than 
20% from their BC03 derived stellar mass, with a me- 
dian CB08 stellar mass that is a factor of ~ 1.5 times 
lower than the BC03 stellar mass. The CB08 models 
also imply ages that are on average 100 Myr younger 
than BC03-derived ages. We will make use of the CB08 
derived ages, stellar masses, and color excesses in the 
subsequent analysis. 

4. DEPENDENCE OF 8 /im LUMINOSITY ON STAR 
FORMATION RATE 

4.1. Computing the Mid-IR and Ha Luminosities 

The 8 /im luminosity, Lg, is computed by fc-correcting 
the 24 fim flux using the average mid-IR spectral shape 
of local starburst galaxies, as discussed in R06. All 90 
galaxies have spectroscopic redshifts and therefore are im- 
mune to the additional 0.3 dex uncertainty introduced by 
photometric redshift errors typical of z ~ 2 star-forming 
galaxies, owing to the complicated mid-IR spectral shape 
and the redshift-sensitivity of the fc-correction (R06). 

Ha luminosities are corrected for extinction based on 
the E(B — V) color excess as derived from the SED fit- 
ting (§ 3). The average color excess is (E(B — V)) — 



0.15 ± 0.10, implying an average UV attenuation of a 
factor of - 4 - 5 (Reddy & Steidel 2004; Reddy et al. 
2006b). Similar average dust obscurations have been ob- 
tained based on stacked X-ray analyses of LBGs at z ~ 3 
(e.g., Seibert et al. 2002; Nandra et al. 2002; Reddy & 
Steidel 2004). Adopting the conversion E(B — V) s tciiar = 
0.4i?(i3-U) nG buiar advocated by Calzetti et al. (2000) re- 
sults in Ha-inferred SFRs that notably overpredict those 
derived from the X-ray and dust-corrected UV for the 
sample presented here(Erb et al. 2006b). Therefore, the 
same value of E(B — V") is assumed for both the UV 
continuum and nebular emission lines. 10 The Calzetti 
et al. (2000) extinction curve then implies a mean Ha 
extinction correction of 1.71 ± 0.40. The uncertainty in 
the observed Ha luminosity is dominated by random er- 
ror in the slit loss correction factor and is taken to be 
w 40% (§ 2.2). 

4.2. X-ray Stacking Analysis 

Recall that our main goal is to determine how mid- 
IR luminosity scales with SFR by cross-correlating the 
Spitzer data with dust-corrected Ha luminosity. We 
have two independent measures of how well the dust- 
corrected Ha luminosity scales with SFR. For instance, 
Erb et al. (2006b) and Reddy k Steidel (2004) demon- 
strate the agreement between dust-corrected UV, dust- 
corrected Ha, and stacked X-ray estimates of the SFRs of 
the same types of galaxies at z ~ 2 that are examined in 
this study. To provide better constraints, and to take ad- 
vantage of the new spectroscopy done after these initial 
studies, we revisited the use of X-ray emission as a proxy 
for SFR in our larger UV spectroscopic sample. To ac- 
complish this, we made several stacks of the deep Chan- 
dra 2 Ms X-ray data in the GOODS-North field (Alexan- 
der et al. 2003) for galaxies in both the Ha and larger 
UV samples, restricting the analysis to those galaxies 
with spectroscopic redshifts. The X-ray stacking analy- 
sis is performed in a manner identical to that specified in 
Reddy & Steidel (2004). The errors on the Ha and UV 
luminosities and SFRs reflect the object-to-object disper- 
sion in Ha/UV luminosity among the stacked galaxies. 
Examination of both the mean and median stacks of the 
X-ray data suggests that there are no galaxies that dom- 
inate the X-ray signal and bias the mean stacked results. 
The results of the stacking analysis are summarized in 
Table 1, along with conversions to SFRs listed in Table 2. 
We will refer to these results throughout the paper. For 

10 Forster Schreiber et al. (2009) find evidence in their sample 
of z ~ 2 galaxies targeted with VTT/SINFONI observations that 
the extinction of the nebular regions may be larger than that for 
the stellar continuum, though their sample includes galaxies with 
somewhat larger SFRs and larger stellar masses than those typical 
of galaxies in our sample. It is possible that the nebular redden- 
ing may on average be larger for galaxies that are forming stars 
at a higher rate and/or hosting an older (less attenuated) stellar 
population (e.g., Epinat et al. 2009). While we note that some 
small amount of additional nebular extinction may be allowable by 
our data after correcting for slit losses (given the f» 40% random 
scatter in such losses), this systematic difference in attenuation is 
likely to be small relative to the scatter in the slit loss correction 
and the scatter in correlation between dust-corrected Ha and dust- 
corrected UV SFRs of ta 0.3 dex. Given that the galaxies (with 
Ha observations) analyzed here are drawn from the same sample 
presented in Erb et al. (2006b), and the fact that these authors 
found no significant evidence for a larger obscuration towards the 
nebular regions in these galaxies, we proceed under the assumption 

that E(B — V)stellar = E(B — U)nebular- 
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TABLE 1 

Stacked X-ray and Ha and UV Luminosities 



Sample 


N a 


-^0. 5 — 2.0 kcV b 


-^2.0 — 10.0 keV C 


J^obs d 
Ha 


j^cor e 
Hq 


J^obs f 
UV 


j^cor g 
UV 


(!) 


( 2 ) 


(3) 






(6) 




(o) 


Ha (All) 


18 


6.6 ± 2.4 


3.4 ± 1.3 


4.3 ± 2.6 


7.4 ± 4.4 


11.5 ± 6.3 


54.4 ± 30.0 


TT /TT 1 i \ h 

Ha (Undct) 


8 


3.3 ± 1.5 


1.8 ± 0.8 


3.5 ± 2.2 


4.9 ± 3.1 


12.2 ± 5.9 


32.4 ± 15.8 


UV (All) 


158 


5.1 ± 1.8 


1.8 ±0.6 






7.6 ±4.1 


41.0 ±22.3 


UV (Undct) h 


116 


2.8 ±1.7 


1.0 ±0.6 






7.4 ±3.9 


29.3 ± 15.4 


UV (0 < SFR < 20) 1 


58 


2.4 ± 1.5 


0.9 ±0.6 






5.9 ±2.6 


15.7 ±7.0 


UV (20 < SFR < 40) 1 


43 


9.6 ±4.6 


3.3 ±1.6 






7.9 ±3.3 


42.9 ± 17.9 


UV (40 < SFR < 80) 1 


32 


12.2 ±5.0 


4.6 ±1.9 






8.4 ±5.3 


81.4 ±31.1 


UV (0 < SFR < 20; Undct) h ' i 


56 


2.0 ± 1.7 


0.7 ±0.6 






5.9 ± 2.7 


15.7 ±7.1 


UV (20 < SFR < 40; Undct) h ' ; 


28 


5.8 ±3.1 


2.3 ±1.2 






8.7 ±3.5 


46.9 ± 18.7 


UV (40 < SFR < 80; Undet) h ' ; 


16 


8.4 ±4.4 


3.9 ±2.1 






8.9 ± 5.6 


78.2 ± 29.5 



a Number of galaxies in sample. 

b Average soft-band X-ray flux, in units of 10 -18 ergs s _1 cm -2 . 

c Average X-ray luminosity in units of 10 41 ergs s , assuming an X-ray photon index of T = 2.0 and Galactic absorption 

column density of 1.6 X 10 20 cm~ 2 (Stark ct al. 1992). 

d Observed Ha luminosity in units of 10 42 ergs s — 1 , uncorrected for dust. 

c Dust-corrected Ha luminosity in units of 10 42 ergs s — 1 . 

f Observed UV luminosity in units of 10 28 ergs s~ 1 Hz~ 1 , uncorrected for dust. 
s Dust-corrected UV luminosity in units of 10 28 ergs s _1 Hz -1 . 

Stack of galaxies not detected at 24 /im. 
1 Stack of galaxies with dust-corrected UV-inferred SFRs in Mq yr" 1 as specified, assuming a Chabrier IMF. 



TABLE 2 

Star Formation Rates and Attenuation from X-ray, Ha, and UV 



Sample a SFR X b SFR°^ b SFRg^ b v) c A^/ ay d SFR°^ b SFR^ b A^ v) c A* v ra; 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



Ha (All) 


37 ±14 19 ±11 


32 ± 19 


~ 1.7 


~ 1.9 


9±5 


42 ± 23 


~4.7 


~ 4.1 


Ha (Undet) 


20 ±9 15 ±10 


22 ± 13 


~ 1.5 


~ 1.3 


9 ± 5 


25 ±12 


~ 2.8 


~ 2.2 


UV (All) 


20 ±7 








6±3 


32 ± 17 


~ 5.3 


~ 3.3 


UV (Undet) 


12±7 








6 ± 3 


23 ± 12 


~ 3.8 


~ 2.0 


UV (0 < SFR < 20) 


10 ±6 








5 ± 2 


12 ±5 


~ 2.4 


~ 2.0 


UV (20 < SFR < 40) 


37 ±18 








6 ± 3 


33 ± 14 


~ 5.5 


~ 6.2 


UV (40 < SFR < 80) 


51 ±21 








7±4 


63 ± 24 


~ 9.0 


~ 7.3 


UV (0 < SFR < 20; Undet) 


8±7 








5 ± 2 


12 ±6 


~ 2.4 


~ 1.6 


UV (20 < SFR < 40; Undet) 


25 ±13 








7±3 


36 ± 15 


~ 5.1 


~ 3.6 


UV (40 < SFR < 80; Undet) 


44 ±23 








7±4 


61 ± 23 


~ 8.7 


~ 6.3 



a Samples are the same as in Table 1. 

b In Mq yr^ 1 , assuming a Chabrier IMF, based on the conversion relations of Ranalli et al. (2003) and Kcnnicutt (1998). 

c Attenuation factor corresponding to E(B — V) assuming the Calzetti et al. (2000) relation. 

d Attenuation factor from the ratio of the X-ray-inferred SFR to the observed Ha and UV SFRs. 



the time being, we point out that the SFRs inferred from 
X-ray emission are broadly consistent with those derived 
from both the UV and Ha emission once the latter two 
are corrected for dust assuming the Calzetti extinction 
curve. Consequently, the mean dust correction derived 
using the Calzetti prescription is similar to the mean 
dust correction obtained by comparing the X-ray SFRs 
with those derived from the observed UV and Ha SFRs, 
even when restricting the sample to those galaxies within 
a specific range of SFR (Figure 3). We will revisit the 
correlation between UV slope and dust attenuation in 
§ 5. Based on the agreement between X-ray and dust- 
corrected Ha and UV inferences of the star formation 



rates, we adopt the dust-corrected Ha luminosity as a 
proxy for the bolometric star formation rate. 

4.3. Correlation Between Ha and 8 fim Luminosity 

Figure 4 shows Lg versus dust-corrected Ha luminos- 
ity for the 90 galaxies in our sample. 11 Focusing on the 

11 For simplicity, wc did not consider the small fraction (si 15%) 
of objects at z ~ 2.0 - 2.6 targeted with NIRSPEC that did not 
yield Ha detections. Roughly half of the galaxies with Ha non- 
detections are also undetected at 24 fim. Including these in our 
analysis does little to alter the correlation between Lg and dust- 
corrected Ha luminosity. Most of the remaining Ha non-detected 
sources had inaccurate astromctry which may have contributed to 
the lack of detection. 




1 10 100 1000 



Fig. 3. — Dust-corrected Ha SFR versus UV SFR (blue symbols; 
data from Erb et al. 2006b) compared to X-ray determined SFRs 
(red symbols). 

61 detected galaxies, we find a probability P ~ 0.00001 
that the Lg and Ha luminosities are uncorrelated, im- 
plying a 4.4 a significance of the correlation. A linear 
fit to the data yields a slope consistent with unity and 
rms dispersion of 0.23 dex. In principle, stacking the 
24 /Ltm data for the 24 /xm-undetected galaxies should 
provide additional constraints for them. However, unlike 
the Poisson-limited X-ray data, there are additional dif- 
ficulties with stacking 24 jam data. The first is that the 
MIPS observations are typically background-limited, re- 
quiring larger samples for a reliable stack. The second is 
that owing to the beamsize of the MIPS observations, we 
must ensure that we include only well-isolated galaxies in 
the stack to avoid contamination from nearby galaxies, 
irrespective of whether the contaminants are detected at 
24 /im. Unfortunately, there are not enough well-isolated 
24 /im-undetected galaxies to obtain a reliable estimate 
of their stacked flux. Of the 29 undetected galaxies, 16 
are at least 3 " away from any IRAC sources. Stacking 
the data for these 16 galaxies results in an upper limit 
only to their median 8 /im luminosity (Figure 4) . In spite 
of these difficulties, we can use a survival analysis to take 
advantage of the censored data (upper limits in Lg) for 
the undetected galaxies to make a more robust measure- 
ment of the overall correlation between Ha and mid-IR 
luminosity. 

There are several well-documented statistical tech- 
niques available for performing linear regression in the 
presence of censored data, including the Buckley-James 
(BJ) and expectation-maximization (EM) algorithms, 
both of which are implemented in the ASURV statistical 
software package (Isobe et al. 1986). Survival analysis as- 
sumes that the presence of an upper limit in a parameter 
(e.g., flux) is independent of the true value of the quan- 
tity being correlated (e.g., luminosity). This is usually 
satisfied in flux-limited surveys where the fundamental 
parameter is luminosity. Moreover, in our case, the 8 /im 
luminosities depend on distance in a highly non-linear 
way because of the strong redshift dependence of the k- 
correction. Finally, the upper limits are determined by 
the depths of the 24 /jm images which are independent 
of the luminosity distribution (i.e., we did not integrate 
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Fig. 4. — 8 fim versus dust-corrected Ha luminosity for the 90 
galaxies with both 24 fim observations and either longslit or nar- 
rowband measurements of the Ho flux. The best-fit linear relation 
for the 61 detected objects (circles) is shown by the dashed line. 
The solid line shows the fit for all galaxies, including the 29 unde- 
tected galaxies (downward-pointing arrows). A stack of the 24 fim 
data for 16 well-isolated undetected galaxies yields a 3 a upper limit 
(large arrow), which is measured from the background dispersion 
in the stacked image. The width of the upper limit is determined 
from the dispersion of Ha luminosities of these 16 galaxies. The 
typical error on individual points is indicated in the panel. 

for a longer period of time on fainter objects). Therefore, 
we assume that the 24 /im upper limits are roughly ran- 
domized with respect to the underlying 8 /im luminosity 
distribution. 

In the EM algorithm, an inital guess is made for Lg for 
the 24 /xm-undetected galaxies based on the sample of 
detections and assuming a normal distribution of residu- 
als about the linear fit. These guesses are then combined 
with the detected sample to produce a new linear fit be- 
tween 8 /im and Ha luminosity. This new fit is then 
used to revise the initial guesses for the luminosities of 
the undetected galaxies, and the procedure is repeated 
until convergence. The BJ method works in a similar 
way, but makes no assumption for the distribution of 
points around the regression line, instead employing the 
Kaplan- Meier distribution derived from the sample (e.g., 
the 24 /im detection probability as a function of Ha lu- 
minosity) . Both methods were used to quantify the rela- 
tion between L(Ha) and Lg using information for all 90 
galaxies, and both methods yielded similar results. In- 
cluding both the detections and non-detections, we find a 
probability P ~ 0.00001 that Lg and L(Ha) are uncorre- 
lated, implying a 4.4 a significance of the correlation. A 
linear fit to the data yields an rms dispersion of 0.24 dex: 

log[L 8 /L ] = 

(1.05 ± 0.11) x log [L{Ha) /ergs s" 1 ] - (34.48 ± 4.26) 
for 10 42 < L{Ha) < 4 x 10 43 ergs s _1 (l) 

Within the errors of the fitted parameters, the linear 
regressions with and without incorporating the upper 
limits are identical over the range of luminosity consid- 
ered here. Note that the correlation between Lg and 
dust-corrected L(Ha) is not solely due to any system- 
atic change in the extinction correction with Lg because 
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Fig. 5. — (Left): Comparison between Lm-Lg relations in the optically thick case and in the case where the emission is both optically-thin 
and thick. Ltr is computed by subtracting the observed UV-derived SFR from the dust-corrected He* SFR (see text). (Right): Comparison 
of several published calibrations between Lg and Lir, with most data taken from Caputi et al. (2007), including those of Reddy et al. 
(2006b); Wu et al. (2005); Caputi et al. (2007); Bavouzet et al. (2008); Chary & Elbaz (2001); Rigby et al. (2008). The calibration and rms 
(random) dispersion derived in our work are denoted by the thick blue line and shaded light blue region. We have plotted these relations 
over only the approximate luminosity range for which they are calibrated. 



we still find a significant 4.0 a correlation between Lg 
and observed L{Ha) (uncorrected for extinction). 

Aside from the statistical analysis, it is worth exam- 
ining further how the 24 /um-undetected galaxies affect 
the overall fit between Lg and L{Ha) based on physical 
grounds. We discussed already the limitations in stack- 
ing the 24 /an data for these galaxies. Therefore, we took 
advantage of the X-ray data for 24 /im-undetected galax- 
ies in the GOODS-N field. An X-ray stack points to an 
average SFR of ~ 20 ± 9 M Q yr _1 , very similar to that 
derived from the dust-corrected Ha (Table 2). Examina- 
tion of the mean attenuation factors suggests also that 
while such 24 /im-undetected galaxies are less dusty than 
the detected galaxies — as expected from the fact that 
24 yum emission traces the dust emission in galaxies — 
the former are not completely dust-free given their rela- 
tively large SFRs and mean attenuation of ~ 1.3 — 1.5. 
Hence, the undetected galaxies cannot have true mid-IR 
luminosities that fall substantially and systematically be- 
low the correlation defined in Eq. 1, and therefore they 
must be described adequately by this correlation. 

Combining the linear fit of Eq. 1 with the relation be- 
tween star formation rate and Ha luminosity (Kennicutt 
1998), which assumes solar metallicity and does not in- 
clude the effects of stellar rotation (e.g., Lcithcrcr 2008), 
we obtain the following: 

log[SFR/M Q yr- 1 ] = 

(0.95 ± 0.10) x log[L 8 /i Q ] - (8.52 ± 0.87) 

for 4 x 10 9 < L s < 2 x 10 11 L Q . (2) 

The SFR in Eq. 2 is a total dust-corrected SFR. In order 
to relate the mid-IR luminosity to total (8 — 1000 /im) in- 
frared luminosity, Lir, we cannot simply combine Eq. 2 
with the Kennicutt relation between Ltr and SFR. This 
is due to the fact that the latter is derived under the 
optically-thick limit in which a vast majority of the bolo- 
mctric luminosity is obscured and emerges in the IR. 
While this is true for most of the galaxies in our sam- 
ple, it will not be the case for fainter ones (e.g., such 



as those undetected at 24 /im) where a large fraction 
of the luminosity may be unobscured (§ 6). The dust- 
corrected Ha-inferred SFR is equivalent to the sum of 
the obscured (IR) and unobscured, or observed, (UV) 
SFRs. Therefore, we subtracted the observed UV SFR 
from the dust-corrected Ha SFR to determine the ob- 
scured SFR and IR luminosity. Taking the observed UV 
component of the SFRs into account, we arrive at the 
following relation between L% and Lir: 

log[£i R /L ] = 

(1.37 ± 0.16) x log[i 8 /i Q ] - (3.01 ± 1.34) 

for 4 x 10 9 < L s < 2 x 10 11 L Q . (3) 

The importance of the correction for the unobscured 
SFR is illustrated in Figure 5 where we compare Eq. 3 
with that obtained assuming optically-thick star forma- 
tion. Galaxies with large 8 /zm dust luminosities have a 
smaller fraction of their luminosity escaping in the UV, 
thus Eq. 3 approaches what we would expect in the case 
where all of the star formation is optically-thick. Al- 
ternatively, galaxies with small dust luminosities have a 
substantial fraction of luminosity emergent in the UV. 
In this case, the obscured SFR is lower than what we 
would have guessed in the optically-thick case. The re- 
sult is that the correlation between and Ltr predicts 
a factor of 4 lower Lir than the optically-thick case at 
I 8 ~3x 10 9 L . 

The equations above summarize the relationship be- 
tween mid-IR luminosity, infrared luminosity, and total 
SFR over the ranges specified above, as derived from the 
correlation between 8 /an and Ha luminosity. 

4.4. Comparison with Other High-Redshift and Local 
Correlations 

Figure 5 summarizes several published determinations 
of the relationship between Lg, and Ltr, compared with 
our determination at z ~ 2 (Eq. 3). There is a gen- 
eral agreement to a factor of w 2 between the con- 
versions over the luminosity range where they overlap 



9 



(10 n ^ Lm, < 10 12 L Q ). Excepting our present deter- 
mination and that of Rigby et al. (2008), all the rela- 
tions shown in Figure 5 were calibrated on galaxies at 
z < 0.7. 12 The calibration presented by Caputi et al. 
(2007) is in rough agreement with our relationship within 
the rms dispersion of the latter. The Wu et al. (2005) 
relations are calibrated on radio and Ha data for local 
galaxies (z < 0.2) drawn from the Spitzer First Look 
Survey (FLS; e.g., Frayer et al. 2006). Finally, the R06 
calibration is derived based on the median L 8 to Lir ra- 
tio of local galaxies in the Elbaz et al. (2002) sample. 
Focusing on the relations calibrated with z ~ 2 data, 
Rigby et al. (2008) report on the correlation between 
Lg and Lir based partly on a sample of 4 lensed galax- 
ies with intrinsic Lir between 10 11 and 10 12 L©. Their 
best-fit relation agrees broadly with ours over the lumi- 
nosity range where the two are calibrated. The value of 
our analysis is that it incorporates the largest sample of 
its kind with Ha, UV, X-ray, and 24 (im observations of 
spectroscopically confirmed typical galaxies at z ~ 2. We 
proceed by adopting Eqs. 1, 2, and 3. 

5. RELATIONSHIP BETWEEN REST-FRAME UV SLOPE 
AND DUST ATTENUATION 

While the previously derived calibrations are useful, 
most high-rcdshift galaxies are too faint or lie at too high 
redshift to be detected directly at 24 fim, thus suggest- 
ing the need for some other proxy for dust extinction. 
Meurer et al. (1999) demonstrated a tight correlation be- 
tween dust attenuation and rest-frame UV slope, (3, for 
a sample of nearby starburst galaxies. An advantage of 
using the UV slope as a proxy for dust obscuration is the 
ability to quantify it for galaxies that are up to two orders 
of magnitude fainter in bolometric luminosity, yet are 
several thousand times more numerous, than the dustiest 
ultraluminous galaxies (ULIRGs) at high redshift. Thus, 
it enables us to quantify the contribution of such galaxies 
to the bolometric luminosity density (e.g., Reddy et al. 
2008; Reddy & Steidel 2009). The Meurer relation is 
widely used to recover the dust attenuation and bolo- 
metric luminosities of high-redshift galaxies, however its 
validity had not been tested directly for such galaxies un- 
til recently. Based on a UV-selected sample of spectro- 
scopically confirmed 24 /im detected galaxies and stacked 
X-ray analysis in the GOODS-N field, R06 demonstrated 
that the local relation appears to hold for galaxies with 
bolometric luminosities (sum of the IR and UV luminosi- 
ties) between 10 11 and 10 12 3 L Q at z ~ 2, luminosities 
typical of L* galaxies at these redshifts. Here, we re- 
evaluate the sensitivity of (3 to dust attenuation using 
the full UV-selected sample of 392 galaxies with 24 /j,m 
data and SED modeling. 

In the subsequent discussion the dust attenuation, 
IRX, is parameterized by the ratio of the infrared lu- 
minosity (computed using Eq. 3) to UV luminosity (un- 
corrected for dust). The intrinsic UV spectral slope (3 
is determined from the G — 1Z color as follows. We gen- 
erated power laws in /(A) oc A* 9 for -2.5 < /3 < 1.0 

12 Bavouzet et al. (2008) test the validity of their results for 
higher redshift galaxies by comparing with stacked 24, 70, and 
160 (im measurements (from the 3 MIPS bands) for galaxies at 
z 1.7, and they find a mean Lg to Lir ratio similar to that of 
z < 0.7 galaxies. Both of their fits with and without the stacked 
data are shown in Figure 5. 



with A/3 = 0.01. These are attenuated for the Lya for- 
est opacity assuming the Madau (1995) prescription and 
multiplied by the G and 1Z transmission filters. The UV 
slope is taken to be the one which gives a G — 1Z color 
closest to the observed value. In the present analysis, 
(3 is simply an observed quantity and is immune to the 
uncertainties associated with SED-fitting (e.g., as would 
be the case with a SED-derived value of E(B — V)). The 
error in UV slope is related directly to the error in color, 
and is typically ~ 0.11. 13 

5.1. Results for Typical Star-Forming Galaxies at z ~ 2 

To aid our discussion, we focus first on the correla- 
tion between UV slope and attenuation for typical star- 
forming galaxies at z ~ 2, which we define formally 
as those galaxies with bolometric luminosities, Lboi = 
Lir + Ltjv < 10 12 L©, and ages older than 100 Myr as 
determined from the broadband SED fitting. These lim- 
its arc adopted to reflect the characteristics of the vast 
majority of star-forming galaxies at z ~ 2. Constraints 
on the bolometric luminosity function imply that a typ- 
ical galaxy at z ~ 2 will have a characteristic luminosity 
L* hol ~ 2x 10 11 L Q (Reddy et al. 2008). Further, our SED 
modeling implies a median age of 360 Myr with a disper- 
sion of ~ 810 Myr, where « 13% have ages younger than 
100 Myr. There are 311 galaxies that under this defini- 
tion are classified as "typical," and 109 of these (35%) 
are detected at 24 /im. These detected galaxies exhibit 
UV slopes that are correlated with attenuation (3.8 o 
significance) with a formal scatter of 0.38 dex about a 
linear fit (Figure 6). The normalization of this correla- 
tion will depend of course on the depth of the 24 /im data, 
as we might expect that the detected galaxies have dust 
attenuations that are larger than those for undetected 
galaxies for a fixed (3. For a fairer representation that is 
robust to the mid-IR observational depth limitations, we 
have stacked the X-ray data for galaxies in the GOODS- 
N field in three bins of (3, irrespective of whether they are 
detected at 24 /xm. We convert the mean X-ray deter- 
mined SFR to an IR luminosity after taking into account 
the fraction of light emerging in the UV (see discussion in 
§ 4.3), and we find mean attenuation factors that exhibit 
the same trend with /3 as the 24 /im-detected galaxies 
(Figure 6). These results are in broad agreement with 
those obtained by stacking the 24 /im data in bins of f3. 

More generally, the mid-IR non-detections exhibit a 
(3 distribution that is statistically unlikely to be drawn 
from the same parent distribution as the mid-IR detected 
galaxies. A two-sided K-S test indicates a probability of 
0.01 that the distribution in (3 for 24 /im detections is 
drawn from the same distribution as that for the 24 /im 
undetected galaxies. While both samples include galax- 
ies over approximately the full range of /3, those that are 
undetected at 24 /im have UV slopes that are systemat- 
ically bluer on average by A/3 ~ —0.1 than the detected 
ones. 

The IRX-/3 trend combined with our observations of 
the (3 distributions for the detected versus undetected 

13 The G-band is affected by the Lya forest only for redshifts 
2 > 2.5. Therefore, statistical fluctuations in the forest will not af- 
fect our determination of /3 as most of the galaxies in the samples 
considered here are at z < 2.5. For the same reason, Lya contami- 
nation of the G-band flux is not a concern in our determination of 
0- 
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Fig. 6. — (Left): Dust attenuation, parameterized as the ratio between infrared and UV luminosity, versus rest-frame UV slope, /3, for 
typical galaxies at z ~ 2 (see text). The typical error in the 24 (im detected points is indicated at the lower right. The filled pentagons 
show our inferences from stacked X-ray data, where the error bars reflect the 1 cr dispersion in the relevant quantity for the stacked galaxies, 
which is typically larger than the formal uncertainty in the stacked X-ray flux. Similarly, the upper limit and open pentagons denote results 
from the 24 fim stacks. The arbitrarily normalized red and blue histograms show the j3 distribution for galaxies undetected and detected, 
respectively, at 24 (im. (Right): Same as left panel for bolometrically luminous galaxies (ULIRGs). Upper limits are shown for 24 /im 
undetected galaxies. 
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Fig. 7. — 24 fim detection fraction as a function of UV spectral 
slope for typical galaxies with ages > 100 Myr and young galaxies 
with ages < 100 Myr in the case of ages derived assuming a Calzetti 
extinction curve (open squares and solid circles) . Errors assume 
Poisson statistics. The dashed lines show the same when we assume 
an SMC extinction curve in deriving the ages. The histograms show 
the arbitrarily normalized j3 distributions (in the case of Calzetti- 
derived ages) for the two subsamples, and a two-sided K-S test 
indicates a probability P < 0.01 that the two are drawn from the 
same parent distribution (see text). 

galaxies implies then that the 24 /xm-undetected galax- 
ies are on average less dusty than the detected galax- 
ies, in accord with expectations. Further, we find that 
the 24 /im detection fraction increases steadily with ft 
for typical star- forming galaxies, implying that galaxies 
with redder spectral slopes are on average more infrared 
luminous (Figure 7). In the next section we show that 
the relationship between ft and luminosity fails for the 
most luminous galaxies at z ~ 2. 

The correlation between attenuation and UV-slope for 
typical L* galaxies at z ~ 2 is not necessarily surpris- 
ing considering our expectations given the correlation 
between mid-IR luminosity and the Ha luminosity cor- 
rected for dust based on the UV slope (Figure 4). Per- 
haps more importantly, the IRX-/3 relation at z ~ 2 



is indistinguishable from that established for local UV- 
starburst galaxies (Meurer et al. 1999), the same relation 
that is almost always used to determine dust-corrected 
luminosities at high redshift. 14 This agreement has been 
noted previously for galaxies selected using other tech- 
niques (e.g., the -Bzi4f-selected sample of Daddi et al. 
2007) as well as UV-selected spectroscopic samples of L* 
galaxies at z ~ 2 (R06; Adelberger & Steidel 2000). Fur- 
ther, the star formation rates derived by correcting the 
observed UV luminosity for dust using the Meurer re- 
lation agree well with the star formation rates inferred 
from the SED-fitting (left panel of Figure 8). This result 
is of course not surprising since the Calzetti extinction 
curve (which is essentially coincident with the Meurer 
relation) is used in the SED-fitting to estimate redden- 
ing and SFRs, both of which are driven primarily by the 
UV slope. The practical utility demonstrated here is that 
the UV-slope can be used to recover the dust attenuation 
of typical high-redshift star-forming galaxies to within a 
1 cr scatter of w 0.4 dex. This strengthens our confidence 
in using the Meurer relation to recover dust attenuation 
from the UV SED of L* galaxies in the absence of longer 
wavelength data. 

5.2. Bolometrically- Luminous Galaxies 

Turning now to bolometrically-luminous galaxies with 
Lboi > 10 12 L Q , we find that they have systematically 
larger IRX ratios than we would have inferred from the 
Meurer relation. This bias is roughly 4 — 5 times larger 
than the random error in the IRX ratio and is generally 
attributed to the fact that significant amounts of star for- 
mation are completely obscured in the UV (e.g., Reddy 
et al. 2006b; Papovich et al. 2006; Chapman et al. 2005; 

14 Note that Meurer et al. (1999) use a different definition of 
IRX than used here. They define IRX as the ratio of FIR-to-UV 
flux density, where Z/tr. ~ 1.75 X Lfir based on the calibration 
of Calzetti et al. (2000). We have applied this correction to the 
Meurer relation before comparing to our data. 
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Fig. 8. — (Left): Dust-corrected UV estimate of the star formation rate (assuming the Meurer relation) versus the SED-derived star 
formation rates. All SFRs assume a Chabrier IMF. (Right): Bolometric SFR, defined as the sum of the IR and UV-derived SFRs, versus 
the SED-derived star formation rates. Orange and red points indicate galaxies with bolometric luminosities Lbol > 10 12 Lp 

Goldader et al. 2002). Hence, the reddening deduced 
from the UV SED tends to be lower than that inferred 
from more direct tracers of the obscured star formation. 

A pertinent question to address is whether our conver- 
sion from 24 /im flux to infrared luminosity is applicable 
for such luminous galaxies because, as the reader will re- 
call, our calibration is based primarily on galaxies with 
^boi i; 10 12 Lq. An examination of the empirically- 
derived fits shown in Figure 5 illustrates that, with 
the exception of the Bavouzet relation, our estimate of 
-Ltr for these luminous objects is similar to what we 
would have predicted from the other relations. Adopting 
the Bavouzet prediction would lower the IRX ratio by 
« 0.3 dex, an amount which is not enough to account 
for the mean 1 dex offset of > 10 12 L© galaxies from the 
Meurer relation. Further, as we show in § 6, the relation- 
ship between bolometric luminosity and dust attenuation 
implies that ultraluminous infrared galaxies (ULIRGs) 
with Lbol ^ 10 12 L Q on average will have fainter observed 
UV luminosities, and thus IRX ratios that are boosted, 
relative to LIRGs. Consequently, the bias of galaxies 
with Lbol > 10 12 L Q to have UV slopes that systemat- 
ically underestimate their attenuation is likely to be a 
physical effect, as opposed to one driven by luminosity- 
dependent biases in the calibrations between mid-IR and 
total IR luminosity. 

Deviation of the these ultraluminous galaxies from the 
Meurer relation can be seen more clearly in Figure 9, 
where we compare the direct tracer of bolometric lumi- 
nosity (sum of infrared and UV luminosities) with that 
obtained when we combine the UV luminosity with the 
Meurer prediction for the dust attenuation. Over the 
luminosity range probed by the data, the Meurer rela- 
tion successfully predicts the bolometric luminosity for 
LIRGs, but significantly underpredicts the luminosity of 
ULIRGs. This naturally translates to a disagreement 
between the SED-derived SFRs and those derived from 
summing the IR and UV-based SFRs for ultraluminous 
galaxies at these redshifts, since the SED-derived SFRs 
assume the Calzetti extinction curve (right panel of Fig- 
ure 8). The critical point is that the Meurer relation can 
be used reliably up to Lbol — 10 12 L Q at z ~ 2, and 



Fig. 9. — Comparison between the direct tracer of bolometric 
luminosity (Lbol = ^IR + ^uv) an d that obtained by combining 
Ltjv with the Meurer prediction of the dust attenuation. 

thus is valid for typical star-forming galaxies at these 
redshifts. 

5.3. Young Galaxies with Inferred Star- Formation Ages 
of < 100 Myr 

We now turn to the 49 galaxies with ages < 100 Myr; 
this young sample constitutes w 13% of the sample of 392 
galaxies. Only 4 of the 49 galaxies are detected directly 
at 24 /im (Figure 10). About three fourths of the young 
galaxies undetected at 24 /im have limits in IRX that 
imply that they lie below the Meurer relation. Stacking 
the 24 /im data results in a conservative 5 a limit in IRX 
that suggests that these young galaxies are in general less 
attenuated than their older counterparts at a fixed value 
of (3. This result is further supported by a stack of the 
X-ray data for the 31 young galaxies in the GOODS-N 
field that places a 5 a upper limit in IRX that is still 
0.5 dex lower than the Meurer expectation. 

These observations suggest that an SMC-like (as op- 
posed to a Calzetti) extinction curve may be more ap- 
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Fig. 10. — Same as Figure 6 for galaxies younger than 100 Myr, where ages are determined from stellar population modeling assuming 
the Calzetti (left panel) and SMC (right panel) extinction curves. The open squares and small downward pointing arrows denote galaxies 
detected and undetected, respectively, at 24 /im. The large light blue and yellow downward pointing arrows denote the 5 a limits from 
stacks of the 24 /jm and X-ray data, respectively, for the < 100 Myr galaxies. Also shown are direct measurements based on CO and 
Spitzer IRS observations of the two lensed LBGs, MS1512— cB58 and the Cosmic Eye, and the form of the SMC extinction curve presented 
in Pettini et al. (1998). There are fewer points in the right panel since ages determined with the SMC curve are generally older than those 
derived with the Calzetti curve. Thus, the fraction of galaxies considered "young" (< 100 Myr) is significantly smaller for the SMC case 
relative to the Calzetti case (see text). 



propriate in describing young galaxies in our sample. 
Because the ages are derived from the stellar popula- 
tion modeling (§ 3) and are therefore degenerate with 
respect to the assumed extinction curve (Calzetti), it 
seems prudent to determine how the inferred ages are 
perturbed under the assumption of an SMC extinction 
curve. Adopting an SMC curve will generally yield older 
ages relative to the Calzetti assumption because in the 
SMC case a smaller fraction of the optical minus near- 
infrared color is attributed to dust and a larger fraction 
is attributed to an older stellar population (e.g., Shapley 
et al. 2001). 
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Fig. 11. — UV luminosity distribution of galaxies with ages less 
than and older than 100 Myr. 



The results obtained by adopting an SMC curve are 



summarized in Figure 10. As expected, the fraction of 
galaxies considered "young" (< 100 Myr) is a factor of 5 
times lower (9/392) than that obtained with a Calzetti 
curve. More relevant to our discussion here is that even 
when considering an SMC curve, we note a significant de- 
parture of the youngest galaxies from the locus of points 
that are characteristic of older galaxies in the IRX-/3 
plane that lie on the Meurer relation. At face value, these 
arguments may seem a bit circular given the degeneracy 
between the assumed extinction law and inferred ages. 
However, both IRX and f3 are determined independently 
of the stellar population modeling and Figure 10 demon- 
strates clearly that irrespective of the choice of extinction 
curve, the IRX-/3 properties of the youngest galaxies de- 
viate from those of older, more typical ones at z ~ 2 (c.f., 
Figure 6). 

5.3.1. 24 fim Detection Fraction and /3 Distribution 

The differences between the young and old populations 
are demonstrated more clearly in Figure 7 where we show 
the 24 ftm detection fraction for the young galaxies versus 
older more typical star-forming galaxies. While the de- 
tection fraction for typical star-forming galaxies increases 
steadily with /3, that same fraction for the young galaxies 
stays essentially flat (with no detected objects). The only 
young galaxies detected directly at 24 ^m are those with 
the reddest spectral slopes. Because the redshift and 
UV luminosity distributions of galaxies in the young and 
typical subsamples are not substantially different (Fig- 
ure 11), the trend between detection frequency and UV 
slope implies that the youngest galaxies are less infrared- 
luminous and less dusty on average than typical galaxies 
at a fixed /?. Note further that the distribution of spectral 
slopes for the older and younger subsamples are unlikely 
to have been drawn from the same parent distribution. A 
K-S test implies a probability P ~ 0.005 that the /3 distri- 
bution of young galaxies is drawn from the same parent 
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distribution as that of the older galaxies. Thus, our anal- 
ysis implies that the youngest galaxies are not only less 
infrared-luminous and less dusty on average than older 
galaxies with the same range in /?, but they also tend to 
have redder spectral slopes than older galaxies. The dif- 
ference in mean f3 between galaxies with ages < 100 Myr 
and > 100 Myr is (/3)<ioo Myr ~ (/3}>ioo Myr = 0.26. The 
apparently contradictory results that the young galax- 
ies have both redder spectral slopes and are on average 
less dusty than typical galaxies may be understood if the 
reddening curve of the young galaxies is "steeper" (i.e., 
less dust is required to redden the continuum) in the UV 
than the typically assumed Meurer/Calzetti relations, as 
alluded to above. To rule out any biases that may under- 
lie the offset of these galaxies from the Meurer relation, 
we must scrutinize the validity of the measures of atten- 
uation for such young galaxies, as we discuss next. 

5.3.2. Measures of Lir in Young Galaxies 

There are several independent reasons to suspect that 
the Lir we compute for these young galaxies are not 
biased by the way in which we are measuring them. 
We consider first the possibility that galaxies with ages 
< 100 Myr have systematically lower Lg/LiR ratios than 
their older counterparts. Specifically, it is thought that 
much of the carbon dust in the ISM is ejected from AGB 
stars on timcscalcs of ~ 200 Myr, whereas the larger 
dust grains responsible for the bulk of the IR emission 
are produced on the much shorter timescale of Type II 
supernovae. Local observations appear to support this 
scenario in the sense that PAH metallicity (Zpah) is ob- 
served to be significantly lower relative to the IR dust 
grain metallicity (Zdust) for metal poor galaxies, but that 
the ratio of the metallicities increases rapidly until an 
equilibrium is reached for gas phase metallicities (Z gas ) 
in excess of w 0.3 Z Q , owing to increased carbon ejection 
by AGBs (Galliano et al. 2008). This delayed injection 
of carbon grains results in a depressed ratio of mid-IR to 
total-IR flux for galaxies younger than 200 Myr. How- 
ever, direct measurements of the mid-IR spectral features 
from Spitzer IRS spectra and long wavelength (> 70 (im) 
constraints on the infrared luminosities of two young 
(< 300 Myr) lensed LBGs at z ~ 3 (MS1512-cB58 
and the "Cosmic Eye") imply Lg/LiR ratios that are 
similar to those of local star-forming galaxies and high- 
redshift submillimeter galaxies (Siana et al. 2008, 2009), 
and older LBGs at z ~ 2 - 3 (Reddy et al. 2006a). In 
addition, an independent deduction based on CO obser- 
vations of MS1512-cB58 (Baker et al. 2001) points to an 
Lir that is not significantly different than what we would 
have predicted from our calibration between L 8 and Lir. 
Further, if we restrict the local correlation between Zpah, 
Zdust, and Z gas to the same dynamic range in gas-phase 
metallicity as observed among LBGs at z ~ 2, 0.3 < 
Zgas < Z (Erb et al. 2006a), we would find a PAH to 
dust grain metallicity ratio that is roughly constant. The 
limited dynamic range in gas-phase metallicity probed by 
the current sample, as well as the direct measurements 
of the PAH to IR ratio in lensed LBGs at z ~ 3, suggests 
that our conversion between L§ and Lir (calibrated pri- 
marily on galaxies older than 100 Myr) should be a rea- 
sonable approximation for the younger galaxies as well. 
In this case, the mid-IR emission observed in the young 
galaxies may be due to some combination of silicate and 



amorphous carbon dust produced in Type II supernovae 
(Todini & Ferrara 2001). 

Finally, the X-ray data provide an independent confir- 
mation of our results. The X-ray emission in starburst 
galaxies arises primarily from shock-heated outflowing 
gas (Strickland et al. 2004; Grimes et al. 2005, 2006) and 
high-mass X-ray binaries (Ghosh & White 2001), both 
of which should be sensitive to star formation on rela- 
tively short timescales (~ 70 Myr) that are comparable 
to the dynamical timescale noted above. Consequently, if 
the true infrared luminosities of the young galaxies were 
larger than what we infer from Ls (and if Lir was con- 
sistent with the Meurer prediction) , then the X-ray stack 
of the 49 young galaxies should have yielded a significant 
detection. Yet, we are able to place a firm 5 a upper limit 
on the the X-ray-inferred IRX of log(LiR/L uv ) < 0.35. 

In summary, galaxies with ages < 100 Myr appear to 
have lower attenuation on average than their UV slopes 
would imply from the Meurer relation (note that the 
Meurer and Calzetti relations track each other closely 
for f3 > —2.0). This inference is based on the lim- 
ited dynamic range in metallicity probed by UV-bright 
(71 < 25.5) LBGs, direct measurements of Ls and Lir 
in at least a couple of young lensed LBGs, stacked X- 
ray measurements, and the 24 /xm detection fraction as 
a function of UV slope. Similar conclusions are reached 
by R06 and Siana et al. (2008, 2009). We have shown 
here that this result appears to apply generically to most 
young < 100 Myr galaxies at z ~ 2. The small fraction 
of young galaxies that have red /? and are detected at 
24 /zm points to a scatter in IRX for young galaxies that 
may be larger than for their older counterparts. 15 Di- 
rect detection of the dust emission of a statistical sample 
of young high-redshift galaxies is required to accurately 
constrain their scatter in attenuation. 

We conclude by noting the following. First, only 10 
of the 90 galaxies used to constrain the relation between 
L§ and L(Ha) have (Calzetti-derived) ages < 100 Myr, 
and redetermining the relation excluding these 10 sources 
does little to affect the overall fit between L 8 and L(Ha). 
In other words, our assumption of the Calzetti relation 
in dust correcting the Ha luminosities for these sources 
minimally affects our conversion between Lg and Lir. 
Second, when computing bolomctric SFRs, the overesti- 
mation of attenuation by applying the Meurer/Calzetti 
relations may be partly compensated for by the fact that 
the UV luminosities in young galaxies will underpredict 
the star formation rate. The latter effect is due to the 
fact that for a constant star formation history the ra- 
tio of O and B stars contributing to the UV continuum 
will stabilize only after the main sequence lifetime of B 
stars of ~ 100 Myr. This underscores the need to correct 
for both a different attenuation and a different conversion 
between UV luminosity and SFR when inferring the total 
star formation rates of young galaxies (ages < 100 Myr) 
at high redshift. 

5.4. Summary 

We have used the relationship between L§, SFR, and 
Lir to examine the variation of UV slope with dust atten- 

15 For example, applying the Meurer relation to one young LBG 
at z = 2.83 ("Wcstphal-MM8") results in a bolometric luminosity 
similar to that inferred from its detection at 850 (im (Chapman & 
Casey 2009). 
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Fig. 12. — (Left): Bolometric luminosity versus dust attenuation for the sample of 392 galaxies with spectroscopic redshifts between 
1.5 < z < 2.6. (Right): Same as left panel but color coded by UV slope fJ and restricted to galaxies with ages > 100 Myr. 



uation for typical galaxies, bolometrically-luminous ones, 
and the youngest galaxies at z ~ 2. We find a significant 
correlation between UV slope and attenuation for the 
vast majority of typical star-forming galaxies at z ~ 2. 
In this case, the UV slope can be used to recover dust at- 
tenuation to within a scatter of 0.4 dex. Those galaxies 
with bolometric luminosities in excess of 10 12 L Q have 
dust obscurations that exceed by almost a decade those 
predicted from the UV slope, owing to the larger fraction 
of obscured luminosity in these galaxies and, as we show 
in § 6, a decrease in observed UV luminosity relative to 
galaxies with moderate bolometric luminosities. 

Finally, for the < 13% of our sample (exact fraction de- 
pends on the extinction law used to model the galaxies' 
photometry; see discussion above) that consists of young 
galaxies with ages < 100 Myr, the local correlation over- 
predicts dust attenuations at a given ft. This implies 
that such young galaxies may on average follow an ex- 
tinction curve that deviates from the usually assumed 
Meurer/Calzetti. These results are similar to what we 
found previously based on a smaller sample of galaxies 
in the GOODS-North field (R06). Here, we have ex- 
panded these initial results by exploring in more detail 
some of the systematics that may give rise to such an off- 
set. We find that the undetected young galaxies and the 
two lensed LBGs from the literature exhibit ft and IRX 
that are consistent with an SMC dust extinction curve 
(e.g., Pettini et al. 1998). This behavior may be due to a 
difference in covering fraction of dust, where the younger 
galaxies have larger covering fractions consistent with a 
foreground screen of dust as described by an SMC-like 
extinction curve. Alternatively, we cannot rule out the 
possibility that the dissimilar extinction curve for young 
galaxies may arise from differences in the dust grain size 
distribution. This may be the case if the dust giving rise 
to the mid-IR emission in young galaxies is produced 
primarily in Type II SNe and contrasts in size and com- 
position from dust produced in lower mass stars (e.g., 
Maiolino et al. 2004; Todini & Ferrara 2001; see also dis- 
cussion in Siana et al. 2009). 



6. 



RELATIONSHIP BETWEEN BOLOMETRIC 
LUMINOSITY AND DUST ATTENUATION 



Employing our measures of the dust obscuration and 
bolometric luminosities for typical star-forming galaxies 
at z ~ 2, we proceed to examine the relationship between 
these two quantities, its implication for the correlation 
between UV and bolometric luminosity, and its redshift 
evolution. 



6.1. Functional Form and ft Dependence 

Based on the sample of 392 galaxies with MIPS obser- 
vations, we find a strong correlation between Lboi and 
dust attenuation with a formal scatter of 0.23 dex about 
a linear fit (Figure 12). The fit to the present data is es- 
sentially identical to the linear fit to a smaller sample in 
the GOODS-N field (R06), despite the differences in the 
conversion between L§ and Ltr used in the current ver- 
sus previous (R06) study. This similarity (and the small 
scatter) is due to the fact that Lboi = Lir+£tjv is highly 
correlated with dust attenuation, Ltr/Ltjv, according to 
our definitions. Nonetheless, we arc confident of our cal- 
ibration between bolometric luminosity and dust atten- 
uation given that Lboi estimated from the sum of the 
IR and UV luminosities is consistent with Lboi (or SFR) 
estimated from X-ray data and from the dust-corrected 
Ha and UV luminosities (§ 4.2). Further, R06 showed 
that the relationship between Lboi and Lir/Ltjv for UV- 
selected samples remains valid for star-forming galaxies 
selected on their rest-frame optical colors and/or submil- 
limeter emission. More recently, Huang et al. (2009) find 
that IRAC-selected ULIRGs at z ~ 1.9 lie on the same 
relation as defined by the typically less luminous UV- 
selected galaxies (LIRGs). These results imply that for 
the observationally-accessible area of the Lboi-LiR/Lyv 
plane, the relationship defined by UV-selected galaxies 
is not substantially different or biased with respect to 
that defined by galaxies selected by other means (e.g., 
optical and IR selections). Formally, our best- fit relation 
between Lboi and dust obscuration is 

log[L bo i/L ] = 
(0.69 ± 0.03) log[L m /L uv ] + (10.91 ± 0.04) 

for L uv > 10 10 L . (4) 
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Fig. 13. — Bolometric luminosity (Lbol) versus observed UV luminosity (Ltjv). Small points and arrows indicate galaxies detected and 
undetected, respectively, at 24 /im. The 24 /an stacked results are shown by the open squares and large arrows (purple). The shaded 
region denotes ±1 cr about the mean relation implied by the correlation between Lbol an d dust attenuation (Eq. 4). The hashed region 
indicates the area excluded by the fact that Lbol must be greater than Ltjv- The dashed horizontal and vertical lines delineate the values 
of LJj V and L£ ol at z ~ 2 (Reddy & Steidel 2009; Reddy et al. 2008), and the thicked dashed line shows the extrapolation of the relation 
to UV- faint galaxies (see text). 



Similar relations between bolometric luminosity, or SFR, 
and dust attenuation have been found at low and high 
redshifts (Wang & Heckman 1996; Adelberger & Steidel 
2000; Reddy et al. 2006b; Buat et al. 2007; Burgarella 
et al. 2008; Buat et al. 2009). 

The variation of spectral slope with attenuation and 
Lbol can be seen in the right panel of Figure 12 where 
points are color coded by f3. Galaxies with red (3 > —1.00 
are on average more attenuated and more bolomctri- 
cally luminous than galaxies with bluer /3. Galaxies with 
— 1.50 < f3 < —1.00 span a relatively larger range of 
attenuation and Lbol, reflecting the fact that many bolo- 
metrically luminous galaxies have a similar range in /? 
(but much larger attenuation) than less luminous galax- 
ies (§ 5.2). Galaxies with the bluest j3 span a narrower 
range in attenuation and Lbol, both being on average 
lower for these galaxies. 



6.2. Variation of Observed UV Luminosity with 
Bolometric Luminosity 

6.2.1. Bolometric Luminosities of UV-Faint Galaxies 

For further insight, we have recast the correlation 
between bolometric luminosity and dust attenuation 
(Eq. 4) in terms of the observed UV luminosity (Ltjv), 



log[Lboi/i Q ] - log 



log[L uv /i Q ; 

2Q(log[L b ol/i ]-b)/a _j_ 



(5) 



where a = 0.69 ± 0.03 and b = 10.91 ± 0.04, as shown in 
Figure 13. The relationship between Lbol and -Ltjv im- 
plies that galaxies with faint UV luminosities are less 
attenuated than their UV-bright counterparts. A re- 
cent study based on deep GOODS and UDF ACS finds 
that UV-faint galaxies at z ~ 2.5 have systematically 
bluer UV slopes (/?) relative to UV-bright ones, imply- 
ing that the former may be less dusty (Bouwens et al. 
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Fig. 14. — 24 (im non-detected fraction as a function of unob- 
scured UV luminosity (Ltjv)- Error bars assume Poisson statistics. 

2009). Further, Reddy et al. (2008) and Reddy & Stri- 
de! (2009) argue based on physical reasons that such UV- 
faint galaxies are unlikely to be as attenuated on average 
as their UV-bright counterparts. For example, significant 
amounts of dust in these galaxies, combined with their 
large number densities as inferred from the steep faint- 
end slope of the UV LF, would result in stellar mass 
densities and an infrared background significantly in ex- 
cess of those measured (Reddy & Steidel 2009). Taking 
this luminosity-dependent dust correction into account, 
Reddy et al. (2008) and Reddy & Steidel (2009) demon- 
strate that such UV faint galaxies dominate the bolomet- 
ric luminosity density at z ~ 2 — 3. 

Our present (larger) data set suggests a significant de- 
pendence between bolometric and UV luminosity, as ev- 
idenced by the trend between 24 /im nondetection frac- 
tion and UV luminosity (Figure 14). Specifically, this 
trend implies that galaxies with faint Ltjv are also on 
average less IR-luminous, in turn suggesting that they 
are less bolometrically-luminous, than their UV-bright 
counterparts. The trend between 24 fim nondetection 
fraction and Ltjv is further supported by a stacking anal- 
ysis. Specifically, we stacked the 24 fim data in bins of 
Ltjv> including both detections and nondctcctions, and 
limited the stack to galaxies with Lboi < 10 12 L Q . The 
latter restriction is imposed because we are interested 
in determining whether a trend exists between dustiness 
and UV luminosity for typical galaxies at z ~ 2 (i.e., 
excluding bolometrically-luminous galaxies with values 
of Ltjv similar to those of galaxies with lower Lboi)- 
These stacked results yield average Lboi consistent with 
the prediction from combining the Meurer relation with 
the trend between UV slope and continuum magnitude 
(e.g., Bouwens et al. 2009), as illustrated in Figure 13. 
Thus, our present sample provides the first direct evi- 
dence (independent of the UV slope) for a trend between 
UV luminosity on the one hand, and bolometric lumi- 
nosity and dust obscuration on the other. 

At first glance, these results run counter to our previ- 
ous analysis that indicated no correlation between dusti- 
ness and UV magnitude for UV-bright (1Z < 25.5) galax- 
ies (Reddy et al. 2008). More specifically, in the previous 



analysis of the GOODS-North field, the stacked values of 
Lir/Ltjv as a function of UV magnitude had uncertain- 
ties that were sufficiently large that we could not rule 
out the possibility of no trend between dustiness and 
UV-magnitude (see Figure 11 of Reddy et al. 2008; see 
also discussion in § 6.5). Note that despite this lack of 
trend observed with the smaller sample, we investigated 
in detail the systcmatics introduced by assuming various 
relations between dustiness and UV luminosity extending 
from UV-bright galaxies (for which empircal constraints 
on the dust obscuration factors exist) to UV-faint galax- 
ies (see Reddy et al. 2008; Reddy & Steidel 2009). The 
variation of IR and bolometric luminosity with UV lu- 
minosity becomes more apparent with our larger sample 
(Figures 13,14, § 6.5). 

6.2.2. Saturation of UV Luminosity 

The shallower-than-unity slope of the relationship be- 
tween Lboi and dust obscuration implies that the UV 
luminosity turns over or "saturates," at which point any 
additional star formation will be optically-thick (Fig- 
ure 13). As we discuss below, this saturation of the UV 
luminosity with increasing SFR likely explains why no 
correlation was found previously between Lboi and Ltjv 
(e.g., as can also be seen by the scatter of the points in 
Figure 13) since these spectroscopic samples probe a rel- 
atively narrow range of UV luminosity around the value 
of LJjy. The saturation point can only be ascertained 
from the functional fit to the Lboi-LiR/Lijv relation and 
should become more apparent with future measurements 
of the dust attenuation of galaxies fainter than our spec- 
troscopic limit. In any case, galaxies routinely found 
in Spitzer surveys at high redshift (e.g., ULIRGs) are 
typically fainter in the UV than moderately luminous 
galaxies found in optical surveys (e.g., LIRGs). As noted 
above, such ULIRGs appear to follow the same Lb r 
dustiness relation as UV-selected galaxies. These results 
can be understood if the dust obscuration, Lir/Ltjv, in- 
creases more rapidly than the increase of total UV lumi- 
nosity as the bolometric luminosity increases. This may 
largely explain why ULIRGs have IRX ratios far in excess 
of the values predicted by their UV spectral slopes based 
on the Meurer relation (§ 5, Figure 9). The turnover in 
Ltjv implies that even in the presence of galaxies with 
very large SFRs at z ~ 2, the observed UV luminos- 
ity will never be brighter than a certain value which, at 
z ~ 2, corresponds to Ltjv ~ 10 11 L Q . 

6.2.3. Implications for the Shape of the UV LF 

The connection between UV and bolometric luminos- 
ity has important implications for the shape of the UV 
luminosity function. In particular, the saturation point 
for the observed UV luminosity occurs approximately at 
the values of L{j V and L£ ol at z ~ 2. This is not sur- 
prising because, by definition, L* is the point brighter 
than which the number density of sources decreases ex- 
ponentially. This effect can also be seen by comparing 
our previous determinations of the UV LF, which is well- 
described by a Schechter function, with the bolometric 
LF, which has a more power-law like shape and not as 
well defined of a "knee" (Reddy et al. 2008; Reddy & 
Steidel 2009). The important point is that the satura- 
tion effect is a direct result of the increasing extinction 
with Lboi- Consequently, dust obscuration must be the 
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Fig. 15. — (Left): Bolometric luminosity versus dust attenuation for our z ~ 2 sample compared with that of local galaxies from Bell 
(2003), Brandl et al. (2006), and Huang et al. (2009). Also shown are lines of constant UV luminosity at the value of ^uv at z = 2 and 
z = (solid lines), and the UV spectroscopic limit at z ~ 2 (dashed line). (Right): Same as left panel, in addition to the best-fit relation 
between bolometric luminosity and dust attenuation for -Ltjv > 10 10 Lq at z ~ 2 (thick solid line) and its extension to UV-faint galaxies 
(thick dashed line). Also shown are lines of constant UV luminosity for Ltjv = 10 88 — 10 108 Lq. 



dominant effect in modulating the bright-end of the UV 
LF, and thus dictating the value of L*, at z ~ 2. 

6.3. Redshift Evolution 

Comparison with the relation derived locally (e.g., 
Buat et al. 2006) shows that L* galaxies at z ~ 2 ex- 
hibit Ltr/ 'Ltjv ratios that are a factor of 5 times lower 
than those of local galaxies with a similar Lboi (Reddy 
et al. 2006b, 2008), implying a redshift evolution in the 
extinction per unit star formation rate (or per unit Lboi)- 
R06 suggested that this may reflect cither an evolution 
in the dust-to-gas ratios as galaxies age or a change in 
the average sizes of the infrared and UV emitting regions 
of galaxies with redshift. In the following, we discuss the 
importance of the Ltjv limit in accurately assessing this 
evolution, a consideration of the Lboi-LiR/Luv relation 
for UV-faint galaxies at high redshift, and the implica- 
tions of these results for the evolution of individual galax- 
ies. 

6.3.1. Dependence on Ltjv Limit 

Quantifying properly the evolution in the Lb r 
Lir/ Xtjv relation requires a careful consideration of the 
Ltjv limit used to compare samples at different redshifts 
(see also the discussion in Buat et al. 2009). The strong 
redshift evolution in Ly V (Reddy & Steidel 2009) im- 
plies that a relatively shallow UV luminosity limit will 
exclude the parameter space in the Lboi-LiR/Luv plane 
where most normal galaxies in the local universe lie (Fig- 
ure 15), causing one to infer a milder redshift evolution 
in the extinction per unit SFR. It is clear that an ac- 
curate comparison of these relations for typical galaxies 
requires a Ltjv limit that is sufficiently faint to detect 
these average (Ly V ) galaxies at all redshifts in question. 
A sufficiently faint Ltjv limit reveals that Ly V galaxies 
at z ~ 2 are a factor of « 10 times less attenuated than 
local galaxies of the same Lboi (Figure 15). 

6.3.2. Consideration of UV-faint Galaxies at High Redshift 



Clearly the Ltjv limit is also an important considera- 
tion in our interpretation of the Lboi-Lin/Luv relation 
at z ~ 2. This limit is defined by the magnitude cut 
of our spectroscopic sample (1Z = 25.5) and corresponds 
to Ltjv ~ 10 10 L Q (Figures 11, 15). This limit is faint 
enough to detect Ly V galaxies at z ~ 2, where L{j V is 
determined not from the present sample alone, but from 
a consideration of the UV luminosity function measured 
using data extending 1.5 mag fainter than our spectro- 
scopic limit (Reddy & Steidel 2009). Therefore, the Lb r 
Lir/Ltjv relation derived here is not just valid for the 
sample analyzed here, but valid for Ly V (z = 2) galaxies 
in general. The saturation of Ltjv discussed above can 
also be seen in Figures 15 and 16 where the relation be- 
tween Lboi and Lir/Luv at z ~ 2 crosses the same line 
of constant Ltjv more than once. 

Let us now consider how the Lboi-LiR/Luv relation 
may change for galaxies fainter than the Ltjv ~ 10 10 L Q 
threshold. Our finding of lower bolometric luminosities 
for UV-faint galaxies (§ 6.2) suggests that extending the 
Ltjv ~ 10 10 L limit to fainter luminosities is unlikely 
to reveal a very large population of UV-faint galaxies 
with Lboi similar to those of UV-bright ones. Therefore 
the Lboi-LiR,/Luv relation is unlikely to be significantly 
broader at Lboi ~ 10 11 — 10 12 L than what we have 
measured from the spectroscopic sample. 

We have illustrated this point by using the relation 
between UV luminosity and j3 (Bouwens et al. 2009) to 
infer the dust attenuation and bolometric luminosities 
using the Meurer relation, thus extending our measure- 
ment of the Lboi-LiR/Luv relation to UV-faint galaxies 
(Figure 15). Examined over a wider range in UV lumi- 
nosity, it becomes clear that the Lboi-LiR/Lyv relation 
at z ~ 2 is defined by a sequence of galaxies with increas- 
ing Lboi and Lir/Ltjv with increasing UV luminosity. 

6.3.3. Evolution of Individual Galaxies 

Figure 16 highlights another important issue. While 
the dustiness of galaxies of a given Lboi increases with 
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Fig. 16. — Same as Figure 15 showing lines of constant Ltjv and 
Lir luminosity. 

decreasing redshift, this should not be miscontrued to 
suggest that a typical z ~ 2 galaxy will evolve to become 
dustier but retain the same Lboi- Clustering analyses of 
L* galaxies at z ~ 2 imply that they evolve to reside 
in the bulges of spirals and low mass ellipticals by z ~ 
(Conroy et al. 2008). Therefore, the evolutionary track of 
any single galaxy in the Lboi-Lirt/Luv plane is unlikely 
to be perpendicular to either axis. As galaxies on average 
fade between z ~ 2 and the present day, reflecting the 
global decline of the star formation rate density (e.g., 
Reddy & Steidel 2009; Reddy et al. 2008; Madau et al. 
1996; Giavalisco et al. 1996; Steidel et al. 1999), there 
will be a corresponding decrease in their Ltjv and £ir, 
and hence Lboi- The critical point is that without taking 
into account the offset between the z ~ 2 and z ~ 
relations, one would wrongly conclude that the dustiness 
observed in local L£ ol galaxies is similar to that observed 
for galaxies of the same bolometric luminosity at z ~ 2, 
when in fact the dustiness is lower at a fixed Lboi at 
higher redshift. 

The correlation between bolometric luminosity and 
dust attenuation simply reflects the relationship between 
star formation rate and gas surface density, i.e., the 
Schmidt law (Zoran et al. 2006). An increase in the 
dust-to-gas ratio with galaxy age would simply shift this 
relationship with redshift, as is observed (Reddy et al. 
2006b; Buat et al. 2007; Zoran et al. 2006). In § 7 we 
revisit the variation of bolometric luminosity with dust 
attenuation in the context of the oxygen abundances of 
starburst galaxies. Geometrical effects are also likely to 
play an important role in the offset between the z ~ 2 and 
z ~ relations. For instance, in Figure 16, many of the 
z ~ galaxies with Lboi ~ 10 10 L Q are bulge-dominated 
spirals with star formation extending on scales of sev- 
eral tens of kpc. If most of the star formation occurs in 
the more metal-poor parts of the galactic disks (Pilyugin 
et al. 2004; Zaritsky et al. 1994; Vila-Costas & Edmunds 
1992), then the measured Z/tr/Ltjv ratios would be sig- 
nificantly lower than what we would infer based on the 
globally-averaged dust-to-gas ratio, including the metal- 
rich component associated with the bulge. 




Fig. 17. — Bolometric luminosity versus dust obscuration of the 
z ~ 2 sample, compared to that of the local IUE sample that was 
used in large part to calibrate the Meurer et al. (1999) relation. 

6.3.4. Comparison with the Local IUE Sample 

The conclusion from the above observations is that 
some combination of dust-to-gas ratio and size evolution 
will drive the observed progression of the LborLrja/Luv 
relation with redshift. A consequence of this evolution 
is that we expect to see an offset in the metallicity- 
luminosity relation with redshift, which we discuss in 
§ 7. In practice, this redshift evolution in the extinc- 
tion per unit SFR implies that rest-UV selection allows 
one to access galaxies with an increasingly larger range in 
L\,o\ with increasing redshift. The increasing UV trans- 
parency with redshift was noted by Adelberger & Steidel 
(2000); Reddy et al. (2006b) for UV-selected samples and 
shown to be true for galaxies selected using other rest- 
optical and submm flux criteria as discussed in Reddy 
et al. (2006b). Daddi et al. (2007) also discuss this is- 
sue in the context of the differences between local and 
high-redshift ULIRGs. The increasing dynamic range in 
bolometric luminosity probed by UV selection is demon- 
strated more directly in Figure 17 where we compare 
typical galaxies at z ~ 2, for which the Meurer relation 
is found to be valid, with many of the same galaxies that 
were used to calibrate the Meurer et al. (1999) relation 
locally; the latter sample is drawn from observations with 
the International Ultraviolet Explorer [IUE] Kinney et al. 
1993; Heckman et al. 1998; Meurer et al. 1999). Galaxies 
that were used to calibrate the Meurer relation locally 
span the same range of IRX (L IR /L UV ) as ^* g a l ax i es 
at z ~ 2, yet the local sample spans luminosities that 
are anywhere from one to two decades less luminous at 
a fixed Lir/Ltjv- The Meurer relation appears to hold 
for typical galaxies at z ~ 2 despite the fact that the 
relation was calibrated on local galaxies that were signif- 
icantly less luminous. This suggests that the UV SED 
alone can be used to recover the total dust attenuation 
in galaxies with progressively larger Lboi at higher red- 
shift, given that, at higher redshift, galaxies are more 
transparent in the UV at a fixed Lboi (Figures 9, 17). 

6.4. Variations in Ltr, and Ltjv &t & Given Lboi 
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The contours of fixed Luv and Lir shown in Fig- 
ure 16 emphasize the relative change in these quantities 
as the dustiness is varied at a fixed iboi at z ~ 2. At 
^boi ^ 2 x 10 11 Lq, stochastic changes in the dust obscu- 
ration (e.g., such as might be expected if the extinction 
is patchy) results in Luv to vary much more than Lir, 
simply because a larger fraction of the total luminosity is 
obscured (e.g., see also Adelberger & Steidel 2000). For 
faint galaxies with L bo i < 10 10 L Q , an adjustment in the 
dust attenuation will do little to alter the UV luminosity 
while Lir changes more dramatically. Of course, none of 
these observations are particularly surprising, given that 
as we adjust Luv there must be a corresponding change 
in _Lir to keep Lboi constant. 

6.5. Comparison with Previous Results 

Several times throughout the discussion above we have 
mentioned the correlation between UV luminosity on the 
one hand, and bolometric luminosity and dust attenua- 
tion on the other. At face value, these results contrast 
with previous investigations that have found no correla- 
tion between UV luminosity and dustiness (Reddy et al. 
2006b; Adelberger & Steidel 2000), as discussed briefly 
in § 6.2.1. Of course, a constant average dust correc- 
tion with UV luminosity naturally leads to lower bolo- 
metric luminosities with decreasing UV luminosity. It is 
clear, however, that if UV-faint galaxies are significantly 
bluer than their brighter counterparts (Reddy et al. 2008; 
Reddy & Steidel 2009; Bouwens et al. 2009), and their 
bolometric luminosities are consistent with the Meurcr 
prediction (stacked points in Figure 13), then their av- 
erage dust obscuration must also be lower (dashed line 
in Figure 15). Hence, their bolometric luminosities will 
be correspondingly lower. We find direct evidence of this 
from the trend between 24 /im detection fraction and UV 
luminosity (Figure 14) , even for the present sample which 
is larger but covers the same dynamic range in Luv as 
previous studies (Reddy et al. 2006b; Adelberger & Stei- 
del 2000). The trend between UV luminosity and UV 
slope was found over a slightly larger dynamic range in 
Luv than is usually represented in spectroscopic samples 
(Bouwens et al. 2009). 

The limited dynamic range and/or smaller samples 
may have contributed to the apparent lack of correla- 
tion between UV luminosity and dustiness observed be- 
fore. There are a couple of other reasons why such a 
correlation may have been difficult to discern in previous 
studies. The first is that the saturation of Luv occurs 
at a value that lies in the range of Luv that is typically 
probed in spectroscopic surveys. In other words, the non- 
monotonic behavior of Luv with bolometric luminosity 
for Luv si 10 11 Lq, folded in with measurement errors, 
washes out the underlying trend between Luv and Lboi- 
The trend between dustiness and UV luminosity becomes 
more apparent once galaxies above the saturation point 
are removed from the analysis (§ 6.2.1). Second, as seen 
in the previous section, the lack of correlation may be due 
to the strong variance of Luv with stochastic changes in 
the dust obscuration for galaxies of moderate (or higher) 
luminosities (e.g., Adelberger & Steidel 2000). It is clear 
that with larger samples, superior (direct) tracers of dust 
emission, and a larger dynamic range, one can obtain a 
more complete picture of how dust obscuration varies 
with other galaxy properties. 



6.6. Summary 

In this section we have considered the functional form 
of the correlation between bolometric luminosity and 
dust obscuration, its consequence for the shape of the 
UV luminosity function at z ~ 2, and its rcdshift evolu- 
tion. For galaxies brighter than our spectroscopic limit 
of Luv ~ 10 10 L , we find a tight positive correlation be- 
tween bolometric luminosity and dust obscuration (Lb r 
Lm/L\jy). Our sample provides the first direct evidence 
for the decrease of infrared luminosity with decreasing 
UV luminosity. Our results suggest that dust obscu- 
ration is likely the dominant effect in modulating the 
bright-end of the UV LF as evidenced by the saturation 
of UV luminosity with increasing star formation rate. We 
demonstrate the importance of the Luv limit in quanti- 
fying the redshift evolution of the LborLm/Luv relation. 
Using a sufficiently faint limit, we find that L* galaxies 
at z ~ 2 have dust obscuration ratios (Lir/Luv) that 
are roughly a factor of 10 times lower than those of local 
galaxies of the same bolometric luminosity. We exam- 
ine this offset in the context of metallicity in the next 
section. 

7. RELATIONSHIP BETWEEN DUST ATTENUATION AND 
METALLICITY ATZ~2 

It is clear from the discussion above that an increase in 
the dust-to-gas ratios with age plays an important role in 
the observed evolution of the correlation between bolo- 
metric luminosity and dust attenuation. More generally, 
because most of the metals in the ISM will be contained 
in (or depleted onto) the same dust grains that attenu- 
ate starlight, we expect that the dust attenuation should 
correlate directly with metallicity and hence stellar mass. 
The near-IR spectroscopic data for LBGs at z ~ 2 have 
been used to constrain their gas-phase metallicities using 
the [Nil] /Ha ratio. Erb et al. (2006a) compute the stellar 
mass-metallicity relation at z ~ 2 in this manner and find 
that the relation is offset from the local one (Tremonti 
et al. 2004) such that at a fixed stellar mass galaxies at 
z ~ 2 arc w 2 times less metal-rich than present-day 
galaxies. This offset likely reflects the larger gas frac- 
tions at a given stellar mass with increasing redshift (see 
discussion in Erb et al. 2006a). 

7.1. Correlation at z ~ 2 

To test our expectation of a correlation between metal- 
licity and dust attenuation, we have investigated the vari- 
ation of Lir/Luv with stellar mass for the sample of 392 
galaxies (Figure 18). Assuming that the 24 /im unde- 
tected galaxies have a similar scatter in Lir/Luv at a 
given mass as the 24 /xm detected galaxies, we find a scat- 
ter between attenuation and stellar mass of ss 0.46 dex 
about a linear fit between the two. The trend between 
attenuation and mass can be visualized more easily by 
examining the mean attenuation in bins of stellar mass 
as inferred from a stacking analysis. To compute the av- 
erage Lir/Luv > we first stacked the 24 //m emission for 
galaxies in each bin of stellar mass to find their aver- 
age Lir luminosity. The average Lir luminosity is then 
combined with the average Luv luminosity for galaxies 
in each bin of stellar mass to compute the Lir/Luv ratio. 
The error in the ratio takes into account the error in the 
stacked estimate of the Lir luminosity and the error in 
the mean of the Luv luminosity. Doing this, we find that 
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Fig. 18. — (Left): Dust attenuation, parameterized as iia/iuVi as a function of stellar mass for 392 galaxies observed with MIPS. 
Galaxies detected and undetected at 24 /im are indicated by the circles and downward pointing arrows, respectively. The average stacked 
values of the attenuation are indicated by the red points. The mass-mctallicity relation from Erb et al. (2006a) is denoted by the orange 
points with the metallity scale indicated on the right-hand axis. The mass-metallicity relation is scaled to match the stacked estimates for 
the attenuation according to Eq. 6. The dashed horizontal line indicates solar metallicity (Asplund et al. 2004). (Right): Measurements 
of attenuation and oxygen abundance (derived from the N2 index; Pettini & Pagel (2004)) for individual star-forming galaxies at z ~ 2. 
Points are color coded to reflect whether they have upper limits in attenuation and/or metallicity. The vertical black error bar denotes 
the random scatter between the [Nil] /Ha ratio and oxygen abundance. The solid red line indicates the relation found by normalizing 
the mass-attenuation relation to match the mass-metallicity relation (left panel). The dashed horizontal line indicates solar metallicity 
(Asplund et al. 2004). 

galaxies with stellar masses > 10 11 M Q are roughly 100 
times more attenuated on average than those with stellar 
masses < 10 9 ' 5 M Q . Because the gas-phase metallicity 
also appears to be significantly correlated with stellar 
mass (and in the same direction as that observed for the 
mass-attenuation relation), it suggests a close connec- 
tion between dust attenuation and gas-phase metallicity, 
in accord with expectations. Relating the mean obscu- 
ration in a given bin of stellar mass with the metallicity 
expected for that stellar mass from the mass-metallicity 
relation, we obtain the following empirical relation be- 
tween obscuration and gas-phase metallicity as measured 
via the oxygen abundance: 



12 + log(0/if ) = 
(0.18 ± 0.03) log(Li R /Luv) + (8.28 ± 0.03) 
for 1 < Lm/iuv < 40. 



(6) 



The upper limit to which the calibration is valid 
(Lir/Ltjv ~ 40) corresponds to the point at which 
[Nil] /Ha begins to saturate for galaxies with close to so- 
lar metallicity (see next section) . This conversion allows 
one to estimate to within 0.54 dex random scatter the 
gas-phase metallicity expected for a given dust attenua- 
tion for L* galaxies at z ~ 2. This scatter includes both 
the scatter in the mass-metallicity relation (Erb et al. 
2006a) and that in the relation between the [Nil] /Ha ra- 
tio and oxygen abundance. This relation can be useful 
in practice because measuring the dust obscuration for 
large numbers of individual galaxies is observationally 
more feasible than measuring gas-phase metallicities; the 
latter requires near-IR spectroscopy at these redshifts. 
Nonetheless, we caution against the use of this relation 
for dust obscurations that are dissimilar to the ones used 
in the calibration, given the very limited dynamic range 



in metallicity probed by the current sample. Figure 18 
shows the comparison between dust attenuation, stellar 
mass, and metallicity. 

Recall that Eq. 6 is derived from the average attenua- 
tion and metallicity of galaxies with a given stellar mass. 
To better gauge how well this relationship holds for indi- 
vidual galaxies, we examined sources in our sample that 
also have published individual metallicity determinations 
(see Table 9 of Law et al. 2009). To remove any potential 
systematics between metallicity calibrations, we consid- 
ered only galaxies where oxygen abudances were inferred 
from the N2 index (Pettini & Pagel 2004). While there 
were only 10 galaxies in our sample that met these re- 
quirements, their dust attenuations and metallicities are 
consistent with our conversion between the two within 
the measurement errors 0.12 dex scatter between 
the individual measurements and the relationship de- 
termined from stacking; Figure 18). Obviously, larger 
samples of galaxies with individual metallicity measure- 
ments, such as those made possible with multi-object 
near-IR spectrographs (e.g., MOSFIRE; McLean et al. 
2008), will enable more detailed investigations of the con- 
nection between dust attenuation and gas-phase metal- 
licity. 

7.2. Comparison with the Local Relations 

We can gain a broader perspective on the z ~ 2 results 
by comparing them to the local trend between metallicity 
and dust obscuration, where locally the galaxy metallic- 
ities are more robustly measured using a variety of cali- 
brations and observations cover a larger dynamic range 
in metallicity. Heckman et al. (1998) find strong positive 
trends of oxygen abundance (O/H) with £ir/ L\jy and 
f3 for a sample of local starburst galaxies (Figure 19). 
They interpret these trends as a reflection of the increas- 
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Fig. 19. — (Left): Oxygen abundance versus UV slope, /3, for galaxies at z ~ 2. For comparison, data for local starburst galaxies from 
Heckman et ai. (1998) are provided, both using the original metallicity calibrations of that study (open circles) and updating their results 
for the [NII]/H«-derived metallicity (filled circles). Their best-fit linear relation (for the original data) is indicated by the short-dashed 
line. The long-dashed line denotes solar metallicity (Asplund et al. 2004). (Right): Same as left panel for oxygen abundance versus dust 
obscuration, Lir/Ltjv- 



ing extinction of UV emission (and reddening of the UV 
slope) as the dust-to-gas ratio increases with gas-phase 
metallicity. 

Before comparing with the z ~ 2 results, we must first 
account for the systematic differences in the way that 
oxygen abundances are inferred locally versus at high 
rcdshift (Kewley & Ellison 2008). As discussed above, 
the abundances for the z ~ 2 sample are derived based 
on the calibration of [Nil] /Ha with O/H (Pettini & Pagel 
2004). Alternatively, the abundances tabulated in Heck- 
man et al. (1998) are based primarily on either direct 
temperature-sensitive (T c ) methods or the R23 calibra- 
tion. For a robust comparison with the z ~ 2 sample, we 
searched the literature to find [Nil] /Ha measurements 
for all of the Heckman et al. (1998) galaxies. The com- 
parison between the metallicities tabulated in Heckman 
et al. (1998) and those derived from the N2 index is 
shown in Figure 19. The two sets of metallicity mea- 
surements agree well up to solar metallicity, at which 
point the [Nil] /Ha indicator saturates due to nitrogen 
becoming the dominant coolant, an effect that has been 
noted elsewhere (Pettini & Pagel 2004; Kewley & Dopita 
2002; Baldwin et al. 1981). 

Having remeasured the metallicities for the local sam- 
ple using the same abundance indicator adopted at high 
rcdshift, we find that the variations in O/H with (3 and 
IRX for L* galaxies at z ~ 2 are in general agreement 
with the local relations (Figure 19). The implication of 
this agreement is that for a fixed metallicity, the z ~ 2 
and z ~ galaxies have similar dust-to-gas ratios, de- 
spite the fact that the high-redshift ones are 10 — 20 
times more bolometrically luminous. It is here that we 
see another manifestation of the redshift evolution of the 
correlation between bolometric luminosity and dustiness. 
As discussed in § 6, z ~ 2 galaxies are significantly more 
luminous at a fixed dust obscuration than local galax- 
ies. We see a similar effect in the context of metallic- 
ity, in the sense that galaxies with a fixed metallicity (or 



dust-to-gas ratio) are significantly more luminous at high 
redshift. This result is a natural expectation given the 
positive correlation between dust attenuation and metal- 
licity (Figure 19). 



7.3. Bolometric Luminosity 



Metallicity Relation at 
2 



The progression of metallicity with increasing obscu- 
ration (Figures 18) and the strong correlation between 
bolometric luminosity and dust attenuation (§ 6) imply 
the existence of a luminosity-metallicity relation at z ~ 2 
(e.g., see also Erb et al. 2006a; Shapley et al. 2004). Fur- 
ther, the similarity in dust-to-gas ratio and difference 
in luminosity suggests a significant redshift evolution in 
the optical luminosity - metallicity relation, as has been 
noted by Erb et al. (2006a); Shapley et al. (2004) in 
the context of z ~ 2 galaxies, and noted elsewhere for 
samples primarily at lower redshift (e.g., Kobulnicky & 
Kewley 2004; Kobulnicky et al. 2003; Maier et al. 2006, 
2004; Salzer et al. 2005; Lara-Lopez et al. 2009). Quan- 
titatively, a LIRG with Lboi ~ 10 11 L Q at z ~ 2 is on 
average a third less metal-enriched than a galaxy with 
the same (LIRG) bolometric luminosity locally. Simi- 
larly, the results suggest that high-redshift galaxies are 
scaled up versions of local starbursts where, for a fixed 
metallicity, the high-redshift ones are significantly more 
luminous than the local ones. 

These observations can be understood if high-redshift 
galaxies are on average less-evolved than local galaxies, 
where the high-redshift ones have larger gas fractions 
and are forming stars in a less metal-rich environment 
(e.g., Erb et al. 2006a; Reddy et al. 2006b). The red- 
shift evolution of both the luminosity-metallicity and the 
luminosity-obscuration relations is then largely reflective 
of the underlying chemical evolution of galaxies as they 
age. As we have seen (§6), this evolution has important 
implications for the relative transparency (or "optical- 
depth" ) of galaxies and their amenability to UV-selection 
at ever-increasing redshifts. 
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8. CONCLUSIONS 

We use ground-based UV imaging, near-IR spec- 
troscopy, and Spitzer MIPS 24 /im imaging of a large 
sample of Lyman Break galaxies to investigate how the 
8 /zm luminosity (Lg) is related to the Ha luminosity 
(L(Ha)), infrared luminosity (Lir), and star formation 
rate of L* galaxies at z <~ 2. Chandra X-ray data are 
used to provide an independent check of the average star 
formation rates derived in this manner. We then use 
the derived relationships to estimate the bolometric lu- 
minosities (Lboi) and dust attenuation of typical galaxies 
at z <~ 2. Our main conclusions are as follows: 

1. Using a sample of 90 Lyman Break Galaxies 
with Ha spectroscopic and narrowband observations and 
MIPS 24 /im imaging, we hnd a tight correlation between 
L$ and L(Ha) with 0.24 dex scatter. We combine this re- 
sult with the Kennicutt relations, taking care to account 
for the unobscured component of the star formation rate, 
to derive relations between L 8 and Lir/SFR. 

2. Based on a larger sample of 392 galaxies with 
MIPS observations, we find that the rest-frame UV slopes 
(/3) of typical star-forming galaxies at z ~ 2 with ages 
> 100 Myr correlate significantly with dust attenuation, 
parameterized by the ratio of infrared-to-UV luminos- 
ity, Lir/Luv- Galaxies with flatter (bluer) f3 are less 
dusty on average than those with redder j3. Further, the 
correlation between /? and dust attenuation is indistin- 
guishable from that established for local UV-starburst 
galaxies (Meurer et al. 1999; Calzetti et al. 2000), the 
latter of which is almost always used to infer the dust 
attenuation of UV-selected galaxies at high redshift. We 
demonstrate here that the local correlation can be used 
to infer the extinction and bolometric luminosities of 
10 10 < L m < 10 12 L© galaxies at z ~ 2 to within a 
scatter of 0.4 dex. Galaxies with the largest bolomet- 
ric luminosities of Lboi ^ 10 12 L Q have bluer j3 than 
their dust attenuations would imply based on the local 
correlation, or the correlation for the vast majority of 
star- forming galaxies at z ~ 2. This effect is likely due 
to the fact the majority of the star formation in these 
bolometrically luminous and heavily attenuated galax- 
ies is optically-thick at UV wavelengths (see also point 5 
below) . 

3. Separately, < 13% of our z ~ 2 sample consists of 
young galaxies with inferred ages of < 100 Myr. Unlike 
their older and more typical counterparts, these young 
galaxies are significantly less attenuated at a given /?, as 
evidenced by their larger non-detection rate at 24 /im 
and their non-detections in the stacked 24 /im and X-ray 
images. These observations suggest that young galax- 
ies may follow an extinction curve that is different than 
the usually assumed Mcurer/Calzctti; the data for the 
young galaxies are consistent with an SMC-like extinc- 
tion curve. If this is the case, then their dust obscuration 
may be up to a factor of 2-3 lower than the values ob- 
tained by assuming the local correlation between j3 and 
dust attenuation. 

4. We verify our previous result that galaxies with 
larger bolometric luminosities are more heavily attenu- 
ated by dust, and galaxies with redder j3 are also more 
bolometrically luminous on average than those with bluer 
p. Comparison with the local correlation between Lboi 
and dust attenuation implies that L* galaxies at z ~ 2 



with iboi ~ 10 11 L© are an order of magnitude less 
dusty than local galaxies with a similar Lboi- Such an 
effect may be related to an increase of dust-to-gas ratio 
as galaxies age. The redshift evolution in the extinction 
per unit SFR implies that in practice UV selection gives 
access to galaxies with an increasingly larger range in 
Lboi at increasing redshift. It also implies that the UV 
SED alone can be used to recover the total dust attenu- 
ation in galaxies with progressively larger Lboi at higher 
redshift, given that, at higher redshift, galaxies are more 
transparent in the UV at a fixed Lboi- 

5. We recast the correlation between bolometric lu- 
minosity and dust attenuation at z ~ 2 in terms of ob- 
served UV luminosity to come to the following conclu- 
sions. First, galaxies with faint UV luminosities are ex- 
pected to be less attenuated than their UV-bright coun- 
terparts. This result is supported by our data that indi- 
cate that UV-faint galaxies are less IR luminous (and less 
bolometrically luminous) than UV-bright ones. Second, 
we observe that galaxies with very large star formation 
rates (e.g., ULIRGs) are less UV luminous per unit SFR 
than galaxies with lower SFRs (e.g., LIRGs), implying 
that the dust covering fraction is likely increasing more 
rapidly with SFR than the intrinsic UV luminosity. This 
effect results in a saturation of the UV luminosity. This 
saturation occurs at the value of L* at these redshifts, 
implying that the bright-end of the UV LF at z ~ 2 is 
likely modulated by dust obscuration. 

6. Motivated by the expectation of a direct correspon- 
dence between extinction and metallicity, we have ex- 
amined the relationship between dust attenuation and 
stellar mass. Galaxies with stellar mass > 10 11 M© are 
almost 100 times more dusty on average than those with 
masses < 10 9 - 5 M Q . The monotonically increasing mass- 
attenuation and mass-metallicity relations imply a close 
connection between attenuation and metallicity, allow- 
ing us to provide an empirical calibration between the 
two. This calibration is verified with a small sample 
of z ~ 2 galaxies with individual metallicity determi- 
nations. Comparison with the local relationship between 
metal abundance and dust obscuration suggests that L* 
galaxies at z ~ 2 have a similar ratio of dust-to-metals 
as local starbursts, despite the high-redshift galaxies' 
bolometric luminosities being a factor of 10 to 20 times 
larger than local galaxies. Our results imply a redshift 
evolution in the luminosity-metallicity and luminosity- 
obscuration relations that reflects the underlying chemi- 
cal evolution of galaxies as they age. 

We have utilized all the data at our disposal, including 
UV imaging, Ha narrowband and spectroscopic observa- 
tions, and Spitzer MIPS 24 /im imaging, to investigate 
the viability of rest-frame 8 /zm luminosity as a star for- 
mation rate and dust indicator. We go on to show how 
these data indicate that the UV slopes of high-redshift 
galaxies are sensitive to dust in the same way as they 
are in the local universe; thus the UV slope is an im- 
portant proxy for inferring the bolometric luminosities 
of high-redshift galaxies in the absence of longer wave- 
length data. Further analysis will benefit from obser- 
vations with near-IR multi-object spectrographs allow- 
ing simultaneous coverage of at least two near-IR bands, 
which will allow for direct measurements of the nebular 
reddening via the B aimer decrement for z ~ 2 galax- 
ies. Future observations with the Herschel Space Tele- 
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scope will no doubt greatly extend our results. These 
new observations will directly sample the longer wave- 
length (rest-frame 30 fim) dust emission of at least lumi- 
nous LIRGs at z ~ 2, and will yield more robust mea- 
surements of the bolometric luminosities of more typical 
galaxies at z ~ 2. 
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